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Caption  Page 

The  WWSSN  station  terms  inferred  from  a  GLM/MLE  joint  inversion  scheme  which  simultane-  14 
ously  inverts  for  the  seismic  source  sizes,  receiver  terms,  as  well  as  the  path  effects.  The  inver¬ 
sion  of  2733  unknown  parameters  is  carried  out  with  13840  signals,  9080  noise  measurements, 
and  1609  dips  from  217  worldwide  explosions  recorded  at  118  selected  WWSSN  stations.  Only 
paths  within  20  and  95  degrees  are  used.  For  each  station,  our  joint  inversion  scheme  regards 
any  component  in  common  across  all  events  as  the  corresponding  “station  term”,  similar  to  the 
Douglas'  (1966)  LSMF  approach.  The  high  correlation  between  the  tectonic  type  and  the  station 
terms  suggests  that  these  empirical  corrections  do  reflect  the  upper  mantle  conditions  under¬ 
neath  the  receivers.  Darkened  stars  represent  some  of  the  nuclear  test  sites  used  in  this  study. 

The  map  showing  the  “pure  propagation  effed"  (top)  and  the  combined  station  amplification  15 
(bottom)  defined  as  the  sum  of  the  receiver  term  (Figure  1)  and  the  path  effect  for  Novaya  Zem- 
lya  explosions.  The  paths  from  Novaya  Zemlya  to  stations  in  North  America  have  systematically 
faster  arrivals  and  smaller  amplitudes,  suggesting  a  profound  defocusing  effect  on  the  first  ar¬ 
rivals;  while  stations  in  Ireland,  Scotland,  Spain,  Bangladesh,  northern  India,  Pakistan,  Korea, 
and  Kenya  report  slow  arrivals  and  large  amplitudes,  suggesting  a  focusing  effect.  Amplitudes 
for  paths  to  Greenland,  Iceland,  Alaska,  Turkey,  Germany,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico, 

Ethiopia,  and  Hawaii,  however,  seem  to  be  controlled  by  the  anelastic  attenuation  with  slow  rays 
also  assodated  with  small  amplitudes,  and  fast  rays  associated  with  large  amplitudes. 

The  map  showing  the  probable  dominating  mechanism  affecting  each  path  from  Novaya  Zemlya.  16 
Positive  symbols  (crosses)  are  those  paths  associated  with  focusing/defocusing,  and  the  nega¬ 
tive  symbols  are  those  dominated  by  the  strong/weak  anelastic  attenuation.  Symbol  size  gives 
a  qualitative  measure  of  the  dominating  mechanism. 

The  crust  beneath  the  Ardic  Ocean  floor  is  composed  of  oceanic  basins  (including  the  Canada  1 7 

Basin,  Makarov  Basin,  Fram  Basin,  and  Nansen  Basin)  which  are  all  located  in  nearly  parallel 
positions  and  are  separated  by  the  Alpha  Ridge,  Lomonosov  Ridge,  and  Arctic  Mid-Ocean  Ridge 
(Nansen-Gakkel  Ridge).  Remarkable  focusing  and  defocusing  effects  of  surface  waves  pro¬ 
pagating  across  the  region  have  been  observed  in  earlier  seismic  studies.  The  P  waves  of  No¬ 
vaya  Zemlya  explosions  to  stations  in  Greenland  and  Iceland  have  to  tunnel  through  the  upper 
mantle  underneath  the  mid-ocean  ridge  axis  (Greenland  Fracture  Zone,  Mohns  Ridge.  Jan  May¬ 
en  Fracture  Zone,  Jan  Mayen  Ridge,  and  Kolbeinsey  Ridge),  where  the  anelastic  attenuation  is 
expected  to  be  very  strong,  and  hence  the  complexity  in  the  P-wave  signal  is  inevitable. 


iv 


List  of  Figures 


Figure  No.  Caption  Page 

5  Scatter  plot  of  3  different  types  of  station  mb  s  for  Novaya  Zemlya  explosion  791018.  The  39  22 

good  recordings,  7  noise,  and  14  clips  are  shown  with  filled  circles,  V-shaped  downward  arrows, 

and  upward  arrows,  respectively.  The  raw  station  mb s  (top)  have  a  standard  deviation  of  0.29 
m.u.  Applying  the  "primary"  station  corrections  reduces  the  scatter  to  0.21  m.u.  Applying  the 
proposed  “secondary”  corrections  to  count  for  the  path  effects  reduces  the  scatter  further  down 
to  0.11  m.u.  The  dashed  lines  around  the  network-averaged  mb  clearly  illustrate  the  remarkable 
reduction  of  fluctuation  across  the  recording  stations.  The  mean  event  mb  itself  is  not 
significantly  changed,  however.  Among  28  Novaya  Zemlya  events  used  in  this  study,  this  event 
has  the  smallest  scatter  in  the  resulting  rr^  g  values.  The  dramatic  reduction  of  variation  from 
m,  to  n\  9  shows  a  factor  of  nearly  2.7,  as  compared  to  the  worst  case  of  about  1 .26  for  the 
event  80101 1  (Figure  9).  Novaya  Zemlya  events  typically  exhibit  a  reduction  factor  around  2. 

6  Scatter  plot  of  3  different  types  of  station  mb  s  for  Novaya  Zemlya  explosion  780927.  Note  that  23 

the  raw  station  mbs  (top)  have  a  standard  deviation  of  0.32  m.u.,  whereas  the  m^g  have  a  stan¬ 
dard  deviation  of  0.11  m.u.  Among  28  Novaya  Zemlya  events  used  in  this  study,  this  event 
shows  the  most  dramatic  reduction  of  variation  with  a  factor  of  2.8. 

7  Same  as  Figure  5  except  for  the  event  870802.  The  variation  reduction  with  a  factor  of  2.2  is  24 

typical  for  Novaya  Zemlya  explosions. 

8  Same  as  Figure  5  except  for  the  event  661027.  Although  this  event  does  not  show  as  dramatic  25 

reduction  in  variation  as  do  other  events,  a  factor  of  1 .4  still  illustrates  the  robustness  of  our  joint 
inversion  scheme.  Note  that  the  path  correction  proposed  in  this  study  not  only  reduces  the  mb 
scatter  at  stations  that  reported  the  good  signals,  but  it  also  improves  the  data  consistency  of 

the  censored  recordings,  as  indicated  by  the  shifting  of  the  clipped  recordings. 

9  Same  as  Figure  5  except  for  the  event  801011  which  is  known  to  be  a  double  explosion.  26 

10  Regressing  the  rr^9(Pa)  on  the  relative  source  size  'F./'F. (671 021)  determined  by  Burger  32 

etai  (1986).  The  uncertainties  in  the  mbs  and  the  (671021)  are  taken  into  account 

through  800  bootstrap  resamplings.  The  darkened  bundle  is  actually  the  collection  of  all  800  re¬ 
gressions,  each  produced  by  a  possible  realization  of  11  perturbed  (mb,  VF„/H'..(671021))  pairs. 

The  95%  confidence  band  (shown  as  2  curves  around  the  darkened  bundle)  is  narrower  near 
the  centroid  and  wider  towards  both  ends,  as  expected.  The  individual  95%  confidence  intervals 
of  the  two  inferred  parameters  (r'.e. ,  the  slope  and  the  intercept  of  the  calibration  curve)  are 
shown  with  the  dashed  line  in  the  scatter  plot  (bottom).  Note  that  the  dashed  rectangle  is  not 
the  joint  90%  confidence  interval,  however,  due  to  the  highly  correlated  nature  of  the  two  param¬ 
eters. 

11  Same  as  Figure  9  except  that  rrhg{Pb)  are  used.  The  Pb  phase  appears  to  fit  the  relative  33 
source  sizes  slightly  better  than  do  Pa  and  P™, . 

12  Same  as  Figure  10  except  the  nVg^max)  are  used.  34 
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SUMMARY 


The  standard  procedure  used  in  estimating  the  source  size  of  underground  nuclear 
explosions  using  mb  measurements  has  been  to  separate  the  station  terms  from  the 
network-averaged  source  terms.  The  station  terms  thus  derived  actually  reflect  the  combina¬ 
tion  of  the  path  effect  and  the  station  effect,  when  only  those  events  in  a  close  proximity  are 
utilized.  If  worldwide  explosions  are  used  in  the  inversion,  then  the  path  effect  tends  to  be 
averaged  out  at  each  station.  In  either  case,  the  effect  due  to  the  propagation  path  alone 
would  not  be  obvious. 

Under  this  research  contract  (F19628-90-C-0158)  with  emphasis  on  the  study  of  seismic 
wave  propagation  in  Eurasia,  we  further  decompose  the  station  amplification  effect  with  a  joint 
inversion  scheme  which  simultaneously  determines  the  seismic  source  size,  the  path  terms, 
and  the  receiver  terms.  Short-period  P- wave  amplitudes  of  217  worldwide  underground 
nuclear  explosions,  including  28  blasts  from  Novaya  Zemlya,  recorded  at  118  WWSSN  sta¬ 
tions  have  been  used  in  one  single  inversion  to  isolate  the  propagation  complexities  affecting 
the  P- wave  amplitudes.  For  all  28  Novaya  Zemlya  events  in  our  WWSSN  database,  the  new 
mb  factoring  procedure  provides  more  stable  mb  measurements  across  the  whole  recording 
network  with  a  reduction  in  the  fluctuational  variation  by  a  factor  of  up  to  3.  Typical  reduction 
factor  in  the  variation  is  about  2.  The  the  worst  case  of  1 .26  turns  out  to  be  a  double  explo¬ 
sion. 


The  inferred  path  terms  are  then  compared  against  the  travel-time  residuals  to  character¬ 
ize  the  propagation  paths.  Our  result  indicates  that  paths  from  the  northern  test  site  in 
Novaya  Zemlya  to  stations  in  North  America  have  systematically  faster  arrivals  and  smaller 
amplitudes,  suggesting  a  profound  defocusing  effect  on  the  first  arrivals;  while  stations  in  Ire¬ 
land,  Scotland,  Spain,  Bangladesh,  northern  India,  Pakistan,  Korea,  and  Kenya  report  slow 
arrivals  and  large  amplitudes,  suggesting  a  focusing  effect.  Amplitudes  for  paths  to  Green¬ 
land,  Iceland,  Alaska,  Turkey,  Germany,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico,  Ethiopia,  and 
Hawaii,  however,  seem  to  be  controlled  by  the  anelastic  attenuation  with  slow  rays  also  asso¬ 
ciated  with  small  amplitudes,  and  fast  rays  associated  with  large  amplitudes.  The  separation 
of  path  effects  from  station  effects  provides  direct  clues  as  how  an  old  fold-belt  structure  (like 
Novaya  Zemlya)  could  "modulate"  the  short-period  P-wave  amplitude  and  travel-time  pat¬ 
terns.  A  strong  correlation  between  the  P-wave  amplitude  and  Lg  detection  at  teleseismic 
distance  is  also  observed. 
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As  a  byproduct  of  this  study,  we  have  derived  the  yield  estimates  of  Novaya  Zemlya 
explosions  based  on  the  path-corrected  mb  values.  Assuming  the  basic  coupling  and  the 
mantle  condition  at  Novaya  Zemlya  are  comparable  to  those  at  Eastern  Kazakhstan,  the  mb 
bias  relative  to  NTS  at  50KT  level  using  the  path-corrected  m^P^)  values  is  inferred  as 
0.25  and  0.36  magnitude  unit  for  Novaya  Zemlya  and  Semipalatinsk,  respectively.  The 
mb{Pmax)  bias  of  0.1 1  between  Semipalatinsk  and  Novaya  Zemlya  could  be  largely  due  to  the 
difference  in  pP  interference  between  these  two  test  sites  (rather  than  the  seismic  coupling). 
The  relative  source  size  determined  by  Burger  etal.  (1986)  and  the  theoretical  ryield  scal¬ 
ing  are  combined  to  extrapolate  our  mb  scaling  to  the  higher  end.  The  resulting  yield  esti¬ 
mates  of  Novaya  Zemlya  explosions  range  from  2  to  2100  KT,  with  peak  values  at  550  KT 
and  65  KT  for  events  before  and  after  1976,  respectively,  which  are  in  reasonable  agreement 
with  those  in  previous  studies. 

Also  included  in  this  report  is  a  complete  listing  of  path-corrected  mb  values  and  yield 
estimates  of  Semipalatinsk  explosions  which  are  used  as  our  baseline  in  calibrating  Novaya 
Zemlya  explosions.  First  motion  of  the  initial  short-period  P  waves  appears  to  be  a  very 
favorable  source  measure  for  explosions  fired  in  hard  rock  sites  underlain  by  the  stable  man¬ 
tle  (such  as  Semipalatinsk). 
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1.  INTRODUCTION 


The  seismic  data  from  the  northern  nuclear  test  site  in  Novaya  Zemlya  (hereby  denoted 
by  NNZ)  are  important  in  part  because  nuclear  tests  at  the  (relatively  more  understood)  Semi- 
palatinsk  test  site  in  Eastern  Kazakhstan  (hereby  denoted  by  KTS)  probably  will  not  be 
resumed,  and  that  several  large  historical  explosions  were  detonated  at  Novaya  Zemlya. 
Many  seismological  problems  associated  with  Novaya  Zemlya  tests  still  need  to  be  studied 
with  various  approaches. 

There  are  fundamental  differences  in  the  energy  flux  in  short-period  P  waves  from 
underground  explosions  at  different  test  sites.  These  differences  reflect  both  near-source 
conditions  (such  as  the  effect  of  local  rock  competence  on  pP  delay  time),  regional  conditions 
(such  as  the  shallow  high-velocity  defo'  jsing  body  beneath  Pahute  Mesa),  and  global  condi¬ 
tions  (such  as  the  deep  mantle  heterogeneity  producing  the  common  pattern  in  the  Nevada 
Test  Site  [NTS]  P  waves  and  P  coda  (Lay  and  Welc,  1987).  Novaya  Zemlya  is  a  very 
interesting  test  site  in  terms  of  the  seismic  propagation  characteristics.  For  instance,  all 
Novaya  Zemlya  explosions  recorded  at  VKA  have  complex  waveforms  (Douglas  etal.,  1973), 
and  similar  (complex)  signals  from  the  same  test  site  have  been  recorded  at  other  stations  in 
North  America  (Butler  and  Ruff,  1980).  We  here  present  our  Novaya  Zemlya  results  based  on 
the  body-wave  magnitude  study  to  explain  some  of  the  complexity  in  propagation  observed 
previously. 

The  body  wave  magnitudes  used  in  developing  and  applying  the  magnitude-yield  rela¬ 
tionship  are  network  mb  values  which  are  some  "average”  of  the  station  mb.  Previously  a 
major  objective  of  magnitude-calculation  research  was  to  determine  the  network  mb  that  is  not 
biased  by  the  sample  truncation  due  to  the  limited  range  of  the  seismometers.  Ringdal  (1978) 
introduced  the  maximum-likelihood  estimator  [MLE]  to  correct  for  the  statistical  bias  inlrcuuced 
by  data  censoring  from  non-detection.  Von  Seggem  and  Rivers  (1978)  pointed  out  the 
importance  of  accounting  for  the  data  censoring  due  to  signal  clipping.  Blandford  and  Shum- 
way  (1982)  derived  the  general  linear  model  [GLM]  in  the  presence  of  data  censoring  using 
the  Expectation  Maximization  [EM]  algorithm.  They  simultaneously  estimated  event  magni¬ 
tudes  and  station  corrections  in  a  maximum-likelihood  sense.  Jih  and  Shumway  (1989)  re¬ 
examined  and  documented  the  GLM  algorithms,  and  they  discussed  the  uncertainty  assess¬ 
ment  in  the  censoring  situation.  It  is  concluded  that  in  the  multi-parameter  linear  regression 
problem  with  censored  data,  the  scaling  of  o/Vdegrees  of  freedom  still  provides  an  extremely 
good  approximation  of  the  uncertainty  associated  with  each  parameter.  In  the  case  of  non 
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censoring,  such  approximation  can  be  proved  to  be  “exact”.  The  methodological  similarities 
and  differences  between  the  iterative  least  squares  [ILSJ  and  the  maximum-likelihood  estima¬ 
tor  [MLE]  were  also  identified  in  Jih  and  Shumway  (1989). 

Jih  and  Wagner  (1991b)  propose  to  compute  the  new  station  magnitude  m29  for  the  i-th 
event  recorded  at  the  j-th  station  as 

m29(i,j)  =  '°9io[A(i.j)/T(i,j)l  +  B(A(i,j))  -  S(j)  -  F(k(i),j) 

where  A(i,j)  is  the  displacement  amplitude  (in  millimicrons)  and  T(i,j)  is  the  period  (in  seconds) 
of  the  P  wave.  The  B(A)  is  the  distance-correction  term.  S(j)  is  the  station  correction,  and 
F(k(i),j)  is  the  path  correction  for  explosions  from  the  k(i)-th  source  region.  The  resulting  new 
magnitude  is  called  m29  to  avoid  confusion  with  the  defined  in  Marshall  etal.  (1979)  that 
corrects  for  the  source-region  attenuation  and  station  terms  solely  based  on  published  Pn 
velocity.  The  path  corrections  determined  in  this  procedure  provide  direct  and  informative 
clues  to  characterize  the  various  propagation  paths.  We  will  also  use  /tt29  extensively 
throughout  this  study  to  characterize  the  scaling  relationship  with  other  magnitude  measure¬ 
ments  and  to  make  the  yield  estimates. 
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2.  PATH-CORRECTED  UNBIASED  NETWORK  mb  ESTIMATOR 

The  conventional  definition  of  station  magnitude  is  computed  as 

mb  =  log10(A/T)  +  B(A)  [1] 

where  A  is  the  displacement  amplitude  (in  nm)  and  T  is  the  predominant  period  (in  sec)  of  the 
P  wave.  The  B(A)  is  the  distance-correction  term  that  compensates  for  the  change  of  P- 
wave  amplitudes  with  distance  (e.g.,  Gutenberg  and  Richter,  1956;  Veith  and  Clawson,  1972). 
mb  in  [1]  is  also  denoted  as  m,  in  Marshall  etal.  (1979).  The  ISC  bulletin  mb  is  just  the  net¬ 
work  average  of  these  raw  station  mb  values  without  any  further  adjustment. 

Consider  NE  explosions  detonated  at  NF  source  regions  that  are  recorded  at  some  or  all 
of  Ns  stations.  The  conventional  GLM  (General  Linear  Model]  or  LSMF  [Least  Squares  Matrix 
Factorization]  network  mb  (Douglas,  1966;  Ericsson,  1971;  von  Seggern,  1973;  Blandford  and 
Shumway,  1982;  Marshall  etal.  ,1984;  Jih  and  Shumway,  1989;  Murphy  etal.,  1989)  is  the 
least-squares  or  maximum-likelihood  network  average  of  the  "station-corrected”  magnitudes: 

m22(i.j)am,(i.j)-S(j)  [2] 

where  S(j)  is  the  "statistical”  or  "empirical"  receiver  correction  at  the  j-th  station.  In  Marshall 
etal.  (1979),  a  priori  information  about  the  P„  velocity  underneath  each  station  is  used  to 
determine  its  associated  “deterministic"  receiver  correction,  S(j),  and  the  resulting  magnitude 
is  called  rr fc.  The  GLM  receiver  corrections,  however,  are  inferred  jointly  from  a  suite  of 
event-station  pairs,  and  no  a  priori  geophysical  or  geological  condition  is  assumed  (and 
hence  the  different  notation  m22).  The  high  correlation  between  the  tectonic  type  and  the 
GLM  station  terms  suggests  that  the  empirical  station  corrections  do  reflect  the  averaged 
upper  mantle  conditions  underneath  the  receivers,  if  the  azimuthal  coverage  at  each  station  is 
broad  enough. 

Jih  and  Wagner  (1991b)  propose  to  compute  the  new  station  magnitude  m29  for  the  i-th 
event  recorded  at  the  j-th  station  as 

m29(i,j)  ■  m,  (i.j)  -  S(j)  -  F(k(i),j)  =  m2  2  (i.j)  -  F(k(i),j)  [3] 

At  the  j-th  station,  F(k(*),j)  is  a  constant  for  all  events  detonated  in  the  same  k-th  "geologically 
and  geophysically  uniform  region".  Partitioning  a  single  nuclear  test  site  into  several 
"regions"  may  be  necessary  in  order  to  account  accurately  for  the  focusing/defocusing  effects. 
This  m28  is  very  similar  to  the  m3  in  Marshall  etal.  (1979)  except  that,  again,  a  priori 
attenuation  information  of  the  source  region  is  used  in  Marshall  et  al.  (1979)  to  determine  the 
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correction  term,  whereas  Jih  and  Wagner  (1991b)  invert  for  the  path  or  near-source  effects 
from  the  data  empirically.  In  other  words,  the  source-region  corrections  proposed  by  Marshall 
et  al.  (1979)  are  constants  (for  all  explosions  in  the  same  source  region)  regardless  of  the 
location  of  the  seismic  stations,  whereas  the  path/near-source  corrections  in  Equation  [3]  are 
highly  dependent  on  the  source-station  paths. 

We  now  examine  briefly  the  fundamental  difference  between  the  present  scheme  (Equa¬ 
tion  [3])  and  the  previous  GLM  schemes.  In  LSMF  and  the  standard  GLM  schemes  ( e.g ., 
Douglas,  1966;  Ericsson,  1971;  von  Seggern,  1973;  Blandford  and  Shumway,  1982;  Marshall 
etal.,  1984;  Lilwall  et  al.,  1988;  Jih  and  Shumway,  1989;  Murphy  etal.,  1989),  it  is  assumed 
that  the  observed  station  mb( i.j)  is  the  sum  of  the  true  source  size  of  the  i-th  event,  E(i),  the 
receiver  term  of  the  j-th  station,  S(j),  and  the  random  noise,  v(i,j): 

m*(i,j)  =  E(i)  +  S(j)  +  v(i,j)  [4] 

The  receiver  term,  S(j),  is  constant  with  respect  to  all  explosions  from  different  test  sites,  and 
hence  it  would  inherently  reflect  the  "averaged"  receiver  effect  —  provided  the  paths  reaching 
the  station  have  broad  azimuthal  coverage.  When  world-wide  explosions  are  used,  the  stan¬ 
dard  deviation  (a)  of  the  noise  v  in  [4]  is  typically  about  0.3  m.u.  or  larger. 

If  LSMF  or  GLM  is  applied  to  events  within  a  smaller  area  of  source  region,  then  the  o  of 
v  in  [4]  could  reduce  to  0.15-0.2  m.u.  However,  the  result  of  such  “single -test-site  GLM" 
approach  should  be  interpreted  or  utilized  cautiously.  The  mb  values  so  determined  are  excel¬ 
lent  estimates  of  “relative  source  size”  for  that  test  site  only.  If  this  “single-test-site  GLM" 
inversion  is  applied  to  several  test  sites  separately,  it  may  not  be  easy  or  obvious  to  find  a 
consistent  baseline  for  “absolute  yield”  estimation  or  immediate  combination  of  the  (inter-site) 
magnitudes,  since  the  recording  network  is  typically  different  from  one  test  site  to  another, 
and  hence  the  station  terms  are  inevitably  inconsistent.  Furthermore,  the  station  terms 
derived  by  the  “single-test-site  GLM"  may  not  necessarily  represent  the  attenuation  under¬ 
neath  the  receiver  side  alone.  They  could  be  contaminated  or  sometimes  even  overwhelmed 
by  the  path/near-source  effects  shared  by  the  explosions  confined  in  a  narrow  azimuthal 
range.  This  could  explain  the  once  puzzling  and  controversial  phenomenon  Butler  and  Ruff 
(1980)  (also  Butler,  1981;  Burdick,  1981)  reported,  namely  that  using  Soviet  explosions  from 
one  test  site  alone  may  fail  to  discern  the  attenuation  differential  between  the  eastern  and 
western  U.S.  There  is  no  doubt,  however,  that  the  GLM  or  LSMF  type  of  methodology  can 
infer  the  station  terms  which  are  strongly  correlated  with  the  upper  mantle  attenuation  under¬ 
neath  the  stations,  provided  the  seismic  sources  have  a  broad  spatial  coverage  as  did  those 
in  North  (1977),  Douglas  and  Marshall  (1983),  Lilwall  and  Neary  (1985),  Ringdal  (1986),  Jih 
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and  Wagner  (1991b),  and  many  others. 

In  the  present  scheme  ([3]),  however,  we  reformulate  the  whole  model  as 

mb(  i.j)  =  E(i)  +  SO)  +  F(k(i),j)  +  v(i.j)  [5] 

where  F<k(i),j)  is  the  correction  term  at  the  j-th  station  for  the  propagation  effect  or  the  near¬ 
source  focusing/defocusing  effect,  which  is  constant  for  all  events  in  the  k-th  “geologically  and 
geophysically  uniform  region".  For  each  seismic  station,  this  F  can  be  regarded  as  its  azimu¬ 
thal  variation  around  the  mean  station  term  S.  However,  as  explained  previously,  it  would  be 
more  appropriate  to  consider  F  the  path  or  near-source  term  because  the  back  azimuths  at 
the  station  could  be  nearly  identical  for  adjacent  test  sites  (such  as  Degelen  and  Murzhik), 
and  yet  the  "F”  terms  could  be  very  different.  By  incorporating  the  F  term  into  the  model,  the 
a  for  world-wide  explosions  is  reduced  to  about  0.2,  roughly  the  same  level  that  which  a 
“single-test-site  GLM"  could  achieve.  Intuitively,  the  present  scheme  (Equations  [3]  or  [5]) 
provides  a  more  detailed  (and  hence  better)  model  than  that  of  Equation  [4]  in  describing  the 
whole  propagation  path  from  the  source  towards  the  receiver.  Simply  put,  Equation  [4]  yields 
a  stronger  fluctuation  in  the  source  terms.  E,  as  well  as  a  larger  standard  deviation  of  v 
because  each  term  in  the  right-hand  side  of  Equation  [4]  would  have  to  “absorb”  part  of  the 
missing  F  term  in  [5].  This  is  exactly  the  same  reason  why  has  smaller  variation  than 

m\  ■ 

Roughly  speaking,  the  procedure  described  in  (3]  has  the  following  advantages: 

•  It  provides  more  stable  mb  measurements  across  the  whole  recording  network,  as  com¬ 
pared  to  the  conventional  GLM  or  LSMF  procedure  which  only  corrects  for  the  station 
terms.  The  reduction  in  the  standard  deviation  of  network  mb  from  m,  to  m^g  could 
reach  a  factor  of  nearly  3  (Figures  1  and  2).  As  a  result,  the  scatter  in  m29  versus 
log(yield)  is  smaller  than  that  for  other  mb,  as  illustrated  by  Jih  and  Wagner  (1991b)  with 
Semipalatinsk  explosions. 

•  The  resulting  network  mb  values  are  not  significantly  different  from  the  GLM  results  (cf. 
Figures  4  through  8).  Thus  if  the  network  mb  values  derived  by  GLM  or  LSMF  are 
unbiased,  so  are  the  refined  results. 

•  The  separation  of  the  path  effect  from  the  station  effect  is  a  crucial  step  to  study  the  vari¬ 
ous  propagation  phenomena,  which  in  turn  would  improve  our  understanding  of  the 
seismic  source  as  well. 

We  have  applied  this  procedure  to  217  worldwide  explosions,  including  28  blasts  from 
northern  subsite  of  Novaya  Zemlya,  and  the  resulting  m2  9  values  are  listed  in  Table  1.  The 
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118  WWSSN  (World  Wide  Standarized  Seismograph  Network]  stations  are  selected  such  that 
each  station  records  10  or  more  good  explosion  signals.  There  are  13840  signals,  9080  noise 
measurements,  and  1609  clips  from  17  test  sites  that  are  used  to  invert  for  the  2733  unknown 
parameters  with  the  maximum-likelihood  approach.  The  standard  deviation  of  v(i,j)  in  [5]  is 
0.196,  as  compared  to  0.294  if  the  conventional  GLM  (Equation  [4])  is  applied  to  the  same 
data  set.  A  complete  list  of  the  event  magnitudes  and  the  path  corrections  for  each  test  site 
can  be  found  in  Jih  and  Wagner  (1992).  Here  we  limit  the  discussion  to  Novaya  Zemlya 
explosions  only. 
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Table  1.  Path-corrected  mb  of  Novaya  Zemlya  Explosions 


Event 

#  of  Signals1 

Magnitudes 

Date 

Site 

Ns  Nn  Nc 

S.E.M.2 

mb(Pa) 

mb(Pb) 

mb(Pmax) 

661027 

NNZ 

56  0  13 

0.024 

6.069 

6.305 

6.451 

671021 

NN  Z 

53  5  3 

0.025 

5.422 

5.610 

5.781 

681107 

NNZ 

58  1  5 

0.024 

5.603 

5.847 

6.031 

691014 

NNZ 

59  2  7 

0.024 

5.768 

5.962 

6.132 

701014 

NNZ 

35  0  22 

0.026 

6.438 

6.631 

6.808 

710927 

NNZ 

23  0  21 

0.029 

6.243 

6.442 

6.579 

720828 

NNZ 

32  0  11 

0.030 

5.983 

6.232 

6.367 

730912 

NNZ 

23  0  21 

0.029 

6.363 

6.679 

6.772 

740829 

NNZ 

25  0  18 

0.030 

6.140 

6.386 

6.569 

751021 

NNZ 

23  0  17 

0.031 

6.104 

6.330 

6.532 

750823 

NNZ 

27  0  12 

0.031 

6.151 

6.389 

6.507 

760929 

NNZ 

27  4  7 

0.032 

5.250 

5.466 

5.600 

761020 

NNZ 

25  34  0 

0.025 

4.258 

4.509 

4.784 

770901 

NNZ 

25  2  2 

0.036 

5.165 

5.461 

5.596 

771009 

NNZ 

18  22  0 

0.031 

4.225 

4.320 

4.537 

780810 

NNZ 

39  3  18 

0.025 

5.412 

5.639 

5.865 

780927 

NNZ 

42  7  10 

0.025 

5.105 

5.362 

5.524 

790924 

NNZ 

39  2  16 

0.026 

5.292 

5.553 

5.739 

791018 

NNZ 

39  7  14 

0.025 

5.305 

5.500 

5.698 

80101 13 

NNZ 

42  4  6 

0.027 

5.188 

5.439 

5.667 

811001 

NNZ 

43  4  5 

0.027 

5.278 

5.512 

5.674 

821011 

NNZ 

32  11  5 

0.028 

5.123 

5.292 

5.444 

830818 

NNZ 

30  4  5 

0.031 

5.316 

5.521 

5.706 

830925 

NNZ 

31  4  5 

0.031 

5.248 

5.467 

5.644 

841025 

NNZ 

22  3  4 

0.036 

5.191 

5.461 

5.624 

870802 

NNZ 

22  3  6 

0.035 

5.318 

5.529 

5.687 

880507 

NNZ 

20  4  1 

0.039 

5.133 

5.311 

5.460 

881204 

NNZ 

15  4  2 

0.043 

5.279 

5.548 

5.677 

1 :  Ns  -  #  of  signals,  Nn  =  #  of  noise  measurements,  Nc  -  #  of  clips. 
2:  standard  error  in  the  mean. 


3:  double  explosion. 
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The  average  offset  between  our  mb(Pb)  values  of  the  11  events  before  1976  and  those 
in  Burger  et  at.  (1986)  (also  measured  from  the  Pb  phase)  is  -0.042±0.017  magnitude  unit 
[m.u.J.  Using  the  RMS  Lg  measured  at  NORSAR  (Ringdal  and  Fyen,  1991;  Ringdal  and 
Marshall,  1989)  as  a  reference,  there  appears  to  be  a  mb(Pmax)-Lg  bias  of  0.1  m.u.  between 
Novaya  Zemlya  and  Eastern  Kazakh.  This  bias  becomes  negligible,  however,  if  m^Pa)  is 
used  to  compare  with  m^Lg)  (instead  of  using  the  conventional  peak-to-peak  amplitudes)  (cf. 
Table  2).  Comparison  with  Nuttli's  (1986b,  1988)  mb(Lg)  gives  very  similar  result  (Table  3). 
The  indifference  of  mb(Pa)-Lg  between  KTS  and  NNZ  indicates  that  perhaps  these  two  test 
sites  have  very  comparable  basic  seismic  coupling;  and  that  the  discrepancy  in  mb(Pmax) 
excitation  (relative  to  Lg)  could  be  due  to  different  pP  interference.  We  will  give  the  yield 
estimates  of  Novaya  Zemlya  explosions  based  on  this  hypothesis  in  a  later  section. 

As  a  side  remark,  there  is  a  significant  difference  in  mb(Pmax)-RMS  Lg  (NORSAR) 
between  southwestern  (SW)  and  northeastern  (NE)  subsites  of  Balapan,  as  discussed  in 
many  other  studies.  The  value  of  0.1 1  m.u.  shown  here  is  slightly  smaller  than  that  of  previ¬ 
ous  studies,  though.  Regressions  with  yields  published  by  Bocharov  etal.  (1989)  show  that 
NE  explosions  have  positive  Lg  residuals  and  negative  mb  residuals;  whereas  SW  explosions 
show  the  opposite  trend  (Jih  and  Wagner,  1991b,  1992).  Thus  it  would  seem  plausible  that 
the  apparent  mb-Lg  bias  could  be  "enhanced"  by  the  negative  correlation  between  mb  and 
Lg  residuals.  It  is  interesting  to  note  the  much  smaller  mb-Lg  bias  when  Pa  is  used  (Table  2). 
A  three-dimensional  geological  model  of  the  Balapan  test  site  by  Leith  and  Unger  (1989) 
shows  a  distinct  difference  between  the  NE  and  SW  portions  of  the  test  site,  with  the  granites 
closer  to  the  surface  and  the  alluvium  thinner  in  the  southwest  (see  also  Bonham  etal., 
1980).  The  thicker  alluvium  layer  in  NE  region  could  increase  the  waveform  complexity  (as 
first  observed  by  Marshall  etal.,  1984)  and  reduce  the  magnitudes  measured  with  Pmax. 
Nevertheless,  the  first  motion  ( i.e .,  Pa)  should  be  least  affected  by  this  factor,  and  therefore  a 
favorable  source  measure  —  so  long  as  it  is  not  contaminated  by  the  microseismic  noise  at 
the  receiver  site. 

Nuttli  (1987,  1988)  suggests  that  there  is  a  mb  bias  of  about  0.2  m.u.  between  Degelen 
and  Balapan,  with  Degelen  explosions  having  even  larger  mb  excitation  (relative  to  Lg).  We 
do  not  see  such  Degelen-Balapan  bias  with  RMS  Lg  measured  at  NORSAR  (Table  2).  The 
Degelen  data  set  alone  is  too  small  for  decisive  conclusion.  However,  if  we  treat  Murzhik  as 
part  of  Degelen,  as  did  Nuttli  (1987),  the  average  mb(Pmax)-PMS  Lg (NORSAR)  bias  between 
Degelen  and  Balapan  is  only  0.02  m.u.,  which  is  insignificant. 
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Table  2.  m2  $-RMS  Lg  (NORSAR)  at  Various  Sites 


Site' 

m b{Pa)-mb(Lg),  # 

mdPb)~mb(Lg)  .  # 

™b(P max)  ,  # 

BSW 

-0.504±0.01 1  11 

-0.22810.011  11 

+0.02310.015  20 

BNE 

-0.565±0.023  8 

-0.30410.014  9 

-0.09210.012  14 

BTZ 

-0.523±0.045  6 

-0.243i0.020  6 

-0.04110.015  14 

Deg 

-0.48410.046  5 

-0.20710.042  5 

+0.01210.034  5 

Mzk 

-0.56210.073  3 

-0.25910.045  3 

-0.04610.032  3 

KTS 

-0.52410.013  33 

-0.25010.010  34 

-0.02610.010  56 

NNZ 

-0.51910.020  14 

-0.29610.023  14 

-0.12110.024  14 

*:  BSW  «=  SW  subsite.  Bala  pan;  BNE  »  NE  subsite,  Balapan;  BTZ  =  transition  zone,  Balapan;  Deg  =  Degelen  Mountain; 
Mzk  =  Murzhik;  KTS  =  all  5  subsites  in  Eastern  Kazakh  combined;  NNZ  =  northern  subsite,  Novaya  Zemlya. 


Table  3.  Tn29-mb[Lg){ Nuttli)  at  Various  Sites 


Site 

™b(Pid-mb(Lg) ,  # 

mdPb)-mb(Lg),# 

mb(Prrax)~mbiLg) ,  # 

BSW 

-0.54410.050  8 

-0.24410.040  8 

-0.03410.028  16 

BNE 

-0.58410.043  6 

-0.30310.039  7 

-0.07510.031  14 

BTZ 

-0.44410.078  5 

-0.18410.057  5 

-0.03310.024  13 

Deg 

-0.52410.124  5 

-0.19610.112  5 

+0.05010.100  5 

KTS 

-0.52910.033  24 

-0.23910.028  25 

-0.03710.01 7  48 

NNZ 

-0.54910.033  24 

-0.32410. 035  24 
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3.  RECEIVER  AND  PATH  EFFECTS  ON  P  WAVES  FROM  NOVAYA  ZEMLYA 


Figure  1  shows  our  receiver  terms  which  are  inferred  jointly  along  with  the  source-size 
estimates  and  path  terms  from  the  worldwide  explosions.  The  receiver  corrections  derived 
with  this  approach  match  the  average  tectonic  structure  underneath  each  station  very  well, 
mainly  due  to  the  broad  coverage  of  azimuths  at  each  station.  Generally  speaking,  the  sta¬ 
tion  terms  are  positive  in  shield  regions  such  as  Australia,  Canada,  India,  and  Scandinavia, 
and  they  are  negative  in  the  east  Africa  rift  valleys,  mid-ocean  ridges  (e.g.,  Iceland  and 
Azores  islands),  island  arcs  (e.g.,  Indonesia,  Japan,  and  Taiwan),  and  Himalaya  Mountain 
Ranges  (Chaman  Fault,  northern  India,  Nepal,  and  Burman  Arc).  Solomon  and  Toksoz 
(1970)  and  many  other  studies  (e.g.,  Evemden  and  Clark,  1970;  Booth  etal.,  1974)  found 
that  for  stations  in  U.S.,  the  attenuation  is  higher  between  the  Rockies  and  Cascades,  and  in 
the  northeastern  U.S.  This  pattern  is  also  observable  in  Figure  1  (see  also  North,  1977).  Jih 
and  Wagner  (1992)  compare  these  empirically  determined  station  terms  with  those  predicted 
by  Marshall  etal.  (1979)  based  on  published  Pn  velocity,  and  find  that  they  match  rather  well. 
As  North  (1 977)  put  it,  it  is  gratifying  that  a  simple  parameter  such  as  mb  can  be  utilized  to 
reveal  the  tectonics.  It  should  be  noted,  however,  that  our  empirical  station  terms  also  include 
the  effect  due  to  the  crustal  amplification  if  such  local  site  effect  is  shared  by  all  ray  paths 
from  different  test  sites  to  a  particular  station.  This  could  be  the  reason  of  a  few  outliers  such 
as  HNR  (Honiara.  Solomon  Islands),  PMG  (Port  Moresby,  East  Papua  New  Guinea),  RAB 
(Rabaul,  New  Britain),  and  BAG  (Baguio  City,  Luzon,  Philippines)  which  do  not  show  negative 
station  terms  as  would  be  expected  from  the  strong  seismicity  in  that  region  (cf.  Figure  1). 
Another  possible  reason  is  that  these  stations  have  relatively  poorer  azimuthal  sampling  in 
our  data  set,  and  hence  the  station  bias  at  these  three  stations  is  not  well  constrained. 

Figure  2  shows  the  map  of  the  "pure  propagation  effect"  (top)  and  the  combined  station 
amplification  (bottom)  defined  as  the  sum  of  the  receiver  term  in  Figure  1  and  the  path  effect 
for  Novaya  Zemlya  explosions.  The  path  term  at  each  station  can  be  regarded  as  the  azimu¬ 
thal  variation  (towards  the  various  source  regions)  relative  to  the  averaged  station 
amplification. 

We  also  applied  the  joint  inversion  scheme  to  the  travel-time  residual  data  set.  Although 
there  are  much  fewer  events  in  our  database  for  which  the  accurate  WWSSN  travel-time  resi¬ 
duals  are  available,  the  multiple  events  from  each  test  site  still  permits  the  reliable  separation 
of  path  and  station  effects.  The  detailed  discussion  of  the  global  result  is  deferred  to  Jih  and 
Wagner  (1992),  and  only  the  result  of  Novaya  Zemlya  is  presented  here.  Figure  3  is  a  map 


showing  the  product  of  the  travel-time  residual  and  the  magnitude  residual  at  each  station. 
Only  the  path  term  of  each  ray  is  used.  Thus  positive  symbols  represent  the  paths  that 
encounter  some  focusing  or  defocusing  structure,  whereas  the  negative  symbols  are  those 
paths  dominated  by  attenuation  mechanism.  The  paths  from  Novaya  Zemlya  to  stations  in 
North  America  have  systematically  faster  arrivals  and  smaller  amplitudes,  suggesting  a  pro¬ 
found  defocusing  effect  on  the  first  arrivals;  while  stations  in  Ireland,  Scotland,  Spain,  Ban¬ 
gladesh,  northern  India,  Pakistan,  Korea,  and  Kenya  report  slow  arrivals  and  large  amplitudes, 
suggesting  a  focusing  effect.  Amplitudes  for  paths  to  Greenland,  Iceland,  Alaska,  Turkey,  Ger¬ 
many,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico,  Ethiopia,  and  Hawaii,  however,  seem  to  be  con¬ 
trolled  by  the  anelastic  attenuation  with  slow  rays  also  associated  with  small  amplitudes,  and 
fast  rays  associated  with  large  amplitudes.  The  crust  beneath  the  Arctic  Ocean  floor  is  com¬ 
posed  of  oceanic  basins  (including  the  Canada  Basin,  Makarov  Basin,  Fram  Basin,  and  Nan¬ 
sen  Basin)  which  are  all  are  located  in  nearly  parallel  positions  and  are  separated  by  the 
Alpha  Ridge,  Lomonosov  Ridge,  and  Arctic  Mid-Ocean  Ridge  (Nansen-Gakkel  Ridge)  (Perry 
et  al. ,  1 986).  Focusing  and  defocusing  effects  of  surface  waves  propagating  across  the  region 
have  been  observed  in  earlier  seismic  studies  (e.g. ,  Zeng  etal.,  1986).  Perhaps  the  upper 
mantle  underneath  these  oceanic  ridges  has  some  complex  features  which  cause  the  strong 
defocusing  effect  across  many  WWSSN  stations  in  North  America  as  observed  in  Figures  2 
and  3. 

Butler  and  Ruff  (1980)  examined  SP  P-wave  amplitudes  of  Novaya  Zemlya  explosions 
recorded  by  WWSSN  stations  in  U.S.  The  lowest  amplitudes  were  found  in  GOL  (Golden, 
Colorado)  and  ALQ  (Albuquerque,  New  Mexico),  with  a  factor  of  4  lower  than  high  amplitudes. 
Our  study  with  enlarged  data  set  recorded  at  a  global  network  show  very  consistent  result. 
The  stations  showing  the  lowest  amplitudes  are  SHK  (Shiraki,  southern  Honshu,  Japan),  KTG 
(Kap  Tobin,  eastern  Greenland),  GOL,  SJG  (San  Juan,  Puerto  Rico  region),  LEM  (Lembang, 
Java),  and  ALQ.  On  the  other  hand,  COP  (Copenhagen,  Denmark),  HLW  (Helwan,  United 
Arab  Republic),  MSH  (Mashhad,  Iran-Turkmenistan  border),  1ST  (Istanbul,  Turkey),  AQU 
(Aquila,  central  Italy),  and  ESK  (Eskdalemuir,  United  Kingdom)  report  the  highest  amplitude 
for  NNZ  explosions  (Table  4).  Note  that  the  station  COP,  which  has  the  largest  combined  sta¬ 
tion  amplification  for  P  wave  from  NNZ,  is  also  used  in  Nuttli's  (1986b,  1988)  mb(Lg)  study  of 
Balapan  and  NNZ  events.  Nuttli  (1988)  obtained  a  smaller  Q  (Q  is  the  elastic  quality  factor) 
for  paths  from  NNZ  to  another  Scandinavian  station  KON  (Konsberg,  southern  Norway)  in  his 
mb(Lg)  study.  It  is  interesting  to  note  that  KON  has  negative  path  term  in  our  analysis  (Table 
4).  The  high  correlation  between  the  signal  strengths  of  teleseismic  P-  wave  recordings  and 
that  of  teleseismic  Lg  detection  is  not  coincidental,  as  similar  phenomena  have  been 
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observed  with  Semipalatinsk  explosions  at  other  stations  as  well  (Jih  and  Wagner,  1992). 


For  explosions  fired  in  fold  belts  such  as  NZ  (in  an  extension  of  the  Urals  fold  belts)  and 
at  Amchitka  site  (in  the  Aleutian  arc),  Douglas  etal.  (1973)  suggest  that  stations  on  great- 
circle  paths  that  lie  along  the  fold  belts  could  show  complex  signals  and  small  amplitudes. 
The  evidence  Douglas  etal.  (1973)  see  with  SI-BC  (Smithers,  British  Columbia)  and  PG-BC 
(Prince  George,  British  Columbia)  recordings  of  the  Amchitka  explosion  LONGSHOT  points  to 
a  low  Q  zone  along  the  fold  belts.  The  WWSSN  station  COL  (College  Outpost,  central  Alaska) 
happens  to  be  located  on  the  great-circle  path  along  both  fold  belts  of  NNZ  and  Aleutian,  and 
it  does  have  anomalously  small  mb  for  Amchitka  events  (Jih  and  Wagner,  1992)  as  well  as 
moderately  small  mb  for  NNZ  events  (Figure  2),  which  appears  to  be  due  to  low  Q  along  the 
propagation  paths  (Figures  3  and  4)  since  the  station  term  itself  is  nearly  zero  (Table  4).  This 
would  seem  to  be  somewhat  similar  to  the  pattern  Jih  and  Wagner  (1991a)  observe  for  the 
Murzhik  explosions  at  azimuths  towards  northwest  and  southeast,  although  a  defocusing 
mechanism  along  the  Chingiz  Fault  near  Murzhik  test  site  would  be  plausible  in  that  case. 

Figures  2,  3  and  Table  4  show  that  WWSSN  stations  in  Greenland  and  Iceland  (GDH, 
Godhavn,  western  Greenland;  KTG,  Kap  Tobin,  eastern  Greenland;  and  AKU,  Akureyri,  Ice¬ 
land)  have  systematically  slow  arrivals  as  well  as  small  amplitudes  for  Novaya  Zemlya  explo¬ 
sions.  which  suggests  that  the  low  Q  mechanism  proposed  by  Douglas  et  at.  (1973)  is  plausi¬ 
ble  for  these  paths.  Figure  4  shows  all  earthquakes  with  mb  >  4.5  within  40°  from  the  North 
Pole.  The  P  waves  from  Matochkin  Shar  of  Novaya  Zemlya  to  these  stations  have  to  tunnel 
through  the  upper  mantle  underneath  the  mid-ocean  ridge  axis  (Greenland  Fracture  Zone, 
Mohns  Ridge,  Jan  Mayen  Fracture  Zone,  Jan  Mayen  Ridge,  and  Kolbeinsey  Ridge),  where 
the  anelastic  attenuation  is  expected  to  be  very  strong.  In  fact,  the  whole  Mohns  Ridge  is  right 
on  the  path  from  Novaya  Zemlya  to  AKU,  and  hence  the  complexity  in  the  P-wave  signal  is 
inevitable.  Chan  et  al.  (1989)  also  observed  extremely  low  Qp  at  WWSSN  station  AKU  with 
surface-wave  recordings.  Nevertheless,  AKU  does  report  simple  waveforms  and  large  ampli¬ 
tudes  for  explosions  fired  in  Balapan,  Degelen,  and  French  Sahara;  for  these  paths  are  not 
affected  by  the  mid-ocean  ridge  (Jih  and  Wagner,  1992). 

Several  studies  have  noted  the  observed  negative  correlation  between  magnitude  and 
complexity:  as  with  earthquake  signals,  the  station  magnitude  of  an  explosions  measured  from 
a  complex  signal  is  smaller  than  that  measured  from  a  simple  signal  (Davies,  1970).  Douglas 
etal.  (1973,  1981)  propose  an  explanation  of  complexity  which  is  similar  to  that  of  Davies 
(1970)  except  that  the  direct  P  (and  pP)  is  assumed  to  be  strongly  attenuated  by  a  region  of 
relatively  low  Q  and  that  the  later  arrivals  that  form  the  coda  of  the  complex  signals  are  in 
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effect  scattered  signals  that  have  traveled  by  the  relatively  high  Q  paths  and  hence  with  little 
attenuation.  On  this  “weak  signal"  explanation  of  complexity,  the  mb  measured  from  a  com¬ 
plex  signal  should  be  less  than  that  computed  from  a  simple  signal  of  the  same  explosion. 
Our  experience  in  measuring  the  WWSSN  film  chips  is  in  general  agreement  with  this  theory. 
Stations  with  negative  combined  path  and  receiver  effects  do  show  complex  waveforms.  LEM 
is  a  notorious  example  for  having  complex  signals  and  small  amplitudes  for  all  Eurasian  test 
sites  we  studied,  including  NNZ.  Nevertheless,  Amchitka  explosions  recorded  at  LEM  are  not 
complex,  despite  that  LEM  is  located  in  a  tectonic  region. 
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Figure  1.  The  WWSSN  station  terms  inferred  from  a  GLM/MLE  joint  inversion  scheme  which  simultane¬ 
ously  inverts  for  the  seismic  source  sizes,  receiver  terms,  as  well  as  the  path  effects.  The  inversion  of 
2733  unknown  parameters  is  carried  out  with  13840  signals,  9080  noise  measurements,  and  1609  clips 
from  217  worldwide  explosions  recorded  at  118  selected  WWSSN  stations.  Only  paths  within  20  and  95 
degrees  are  used.  For  each  station,  our  joint  inversion  scheme  regards  any  component  in  common 
across  all  events  as  the  corresponding  “station  term”,  similar  to  the  Douglas’  (1966)  LSMF  approach 
The  high  correlation  between  the  tectonic  type  and  the  station  terms  suggests  that  these  empirical 
corrections  do  reflect  the  upper  mantle  conditions  underneath  the  receivers.  Darkened  stars  represent 
some  of  the  nuclear  test  sites  used  in  this  study. 
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Figure  2.  The  map  showing  the  "pure  propagation  effect"  (top)  and  the  combined  station  amplification 
(bottom)  defined  as  the  sum  of  the  receiver  term  (Figure  1)  and  the  path  effect  for  Novaya  Zemlya 
explosions.  The  paths  from  Novaya  Zemlya  to  stations  in  North  America  have  systematically  faster 
arrivals  and  smaller  amplitudes,  suggesting  a  profound  defocusing  effect  on  the  first  arrivals;  while  sta¬ 
tions  in  Ireland,  Scotland,  Spain,  Bangladesh,  northern  India,  Pakistan,  Korea,  and  Kenya  report  slow 
arrivals  and  large  amplitudes,  suggesting  a  focusing  effect.  Amplitudes  tor  paths  to  Greenland,  Iceland, 
Alaska,  Turkey,  Germany,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico,  Ethiopia,  and  Hawaii,  however,  seem  to 
be  controlled  by  the  anelastic  attenuation  with  slow  rays  also  associated  with  small  amplitudes,  and  fast 
rays  associated  with  large  amplitudes. 
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Figure  3.  The  map  showing  the  probable  dominating  mechanism  affecting  each  path  from  Novaya 
Zemfya.  Positive  symbols  (crosses)  are  those  paths  associated  with  focusing/defocusing,  and  the  nega¬ 
tive  symbols  are  those  dominated  by  the  strong/weak  anelas'ic  attenuation.  Symbol  size  gives  a  quali¬ 
tative  measure  of  the  dominating  mechanism. 


Figure  4.  The  crust  beneath  the  Arctic  Ocean  floor  is  composed  of  oceanic  basins  (including  the 
Canada  Basin,  Makarov  Basin,  Fram  Basin,  and  Nansen  Basin)  which  are  all  located  in  nearly  parallel 
positions  and  are  separated  by  the  Alpha  Ridge,  Lomonosov  Ridge,  and  Arctic  Mid-Ocean  Ridge 
(NanservGakkel  Ridge).  Remarkable  focusing  and  defocusing  effects  of  surface  waves  propagating 
across  the  region  have  been  observed  in  earlier  seismic  studies.  The  P  waves  of  Novaya  Zemlya 
explosions  to  stations  in  Greenland  and  Iceland  have  to  tunnel  through  the  upper  mantle  underneath  the 
mid-ocean  ridge  axis  (Greenland  Fracture  Zone,  Mohns  Ridge,  Jan  Mayen  Fracture  Zone,  Jan  Mayen 
Ridge,  and  Kolbeinsey  Ridge),  where  the  anelastic  attenuation  is  expected  to  be  very  strong,  and  hence 
the  complexity  in  the  P-  wave  signal  is  inevitable. 
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Table  4.  Receiver  and  Path  Effect  on  mb  of  Novaya  Zemlya  Events 

Station 

Lon 

Lat 

Rev 

Path 

Rcv+Path 

SHK 

132.678 

34.532 

-0.250 

-0.473 

-0.723 

KTG 

-21.983 

70.417 

-0.254 

-0.321 

-0.575 

GOL 

-105.371 

39.700 

-0.237 

-0.312 

-0.549 

SJG 

-66.150 

18.112 

-0.173 

-0.352 

-0.525 

LEM 

107.617 

-6.833 

-0.417 

-0.066 

-0.483 

ALQ 

-106.457 

34.943 

-0.199 

-0.275 

-0.474 

GDH 

-53.533 

69.250 

-0.155 

-0.275 

-0.430 

DUG 

-112.813 

40.195 

0.074 

-0.446 

-0.372 

BHP 

-79.558 

8.961 

-0.056 

-0.305 

-0.361 

CHG 

98.977 

18.790 

-0.234 

-0.125 

-0.359 

QUE 

66.950 

30.188 

-0.484 

0.161 

-0.323 

MAT 

138.207 

36.542 

-0.256 

-0.056 

-0.312 

ANP 

121.517 

25.183 

-0.139 

-0.164 

-0.303 

BEC 

-64.681 

32.379 

-0.120 

-0.176 

-0.296 

MSO 

-113.941 

46.829 

-0.091 

-0.258 

GUA 

144.912 

13.538 

-0.232 

0.006 

-0.226 

KON 

9.598 

59.649 

0.000 

-0.212 

-0.212 

OGD 

-74.596 

41.088 

-0.164 

-0.039 

-0.203 

FVM 

-90.426 

37.984 

0.069 

-0.266 

-0.197 

LON 

-121.810 

46.750 

-0.121 

-0.074 

-0.195 

AKU 

-18.107 

65.687 

-0.053 

-0.108 

-0.161 

WES 

-71.322 

42.385 

-0.228 

0.067 

-0.161 

DAV 

125.575 

7.088 

-0.040 

-0.101 

-0.141 

EPT 

-106.506 

31.772 

-0.070 

-0.069 

-0.139 

POO 

73.850 

18.533 

-0.004 

-0.135 

-0.139 

COR 

-123.303 

44.586 

0.161 

-0.297 

-0.136 

NAT 

-35.033 

-5.117 

0.128 

-0.250 

-0.122 

TUC 

-110.782 

32.310 

-0.051 

-0.067 

-0.118 

CMC 

-115.083 

67.833 

-0.270 

0.162 

-0.108 
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Table  4.  Receiver  and  Path  Effect  on  mb  of  Novaya  Zemlya  Events 

Station 

Lon 

Lat 

Rev 

Path 

Rcv+Path 

KBL 

69.043 

34.541 

-0.188 

0.087 

-0.101 

BKS 

-122.235 

37.877 

0.104 

-0.202 

-0.098 

GSC 

-116.805 

35.302 

0.022 

-0.113 

-0.091 

HKC 

114.172 

22.304 

-0.087 

-0.001 

-0.088 

UNM 

-99.178 

19.329 

-0.060 

-0.028 

-0.088 

PTO 

-8.602 

41.139 

-0.193 

0.118 

-0.075 

COL 

-147.793 

64.900 

-0.002 

-0.068 

-0.070 

LUB 

-101.867 

33.583 

0.135 

-0.204 

-0.069 

SCP 

-77.865 

40.795 

-0.060 

0.001 

-0.059 

LPS 

-89.162 

14.292 

-0.109 

0.053 

-0.056 

SNG 

100.620 

7.173 

0.005 

-0.044 

-0.039 

ATL 

-84.338 

33.433 

0.058 

-0.093 

-0.035 

NIL 

73.252 

33.650 

-0.069 

0.059 

-0.010 

GEO 

-77.067 

38.900 

-0.006 

0.026 

0.020 

NAI 

36.804 

-1.274 

-0.110 

0.141 

0.031 

KOD 

77.467 

10.233 

0.177 

-0.138 

0.039 

STU 

9.195 

48.772 

-0.001 

0.041 

0.040 

SEO 

126.967 

37.567 

-0.125 

0.170 

0.045 

SHL 

91.883 

25.567 

-0.081 

0.126 

0.045 

AAE 

38.766 

9.029 

-0.290 

0.342 

0.052 

JCT 

-99.802 

30.479 

0.095 

-0.041 

0.054 

TRI 

13.764 

45.709 

-0.193 

0.247 

0.054 

SDB 

13.572 

-14.926 

-0.049 

0.109 

0.060 

LOR 

3.851 

47.267 

0.008 

0.081 

0.089 

FLO 

-90.370 

38.802 

-0.100 

0.195 

0.095 

DAL 

-96.784 

32.846 

0.266 

-0.164 

0.102 

MDS 

-89.760 

43.372 

-0.091 

0.202 

0.111 

KIP 

-158.015 

21 .423 

-0.041 

0.139 

EIL 

34.950 

29.550 

-0.067 

0.219 

0.152 
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Table  4.  Receiver  and  Path  Effect  on  mb  of  Novaya  Zemlya  Events 


Station 

Lon 

Lat 

Rev 

Path 

Rcv+Path 

JER 

35.197 

31.772 

-0.014 

0.178 

0.164 

MAL 

-4.41 1 

36.728 

-0.056 

0.221 

0.165 

OXF 

-89.409 

34.512 

0.169 

-0.004 

0.165 

NDI 

77.217 

28.683 

0.100 

0.080 

0.180 

BAG 

120.580 

16.411 

0.157 

0.187 

BLA 

-80.421 

37.211 

0.022 

0.173 

0.195 

VAL 

-10.244 

51.939 

-0.024 

0.227 

0.203 

MAN 

121.077 

14.662 

0.367 

-0.160 

0.207 

CAR 

-66.928 

10.507 

0.209 

0.012 

0.221 

TAB 

46.327 

38.068 

0.290 

-0.067 

0.223 

ATU 

23.717 

37.972 

0.171 

0.063 

0.234 

TRN 

-61.403 

10.649 

0.140 

0.108 

0.248 

PDA 

-25.663 

37.747 

0.043 

0.246 

0.289 

SHI 

52.520 

29.638 

0.120 

0.170 

0.290 

BUL 

28.613 

-20.143 

-0.004 

mum 

0.308 

RCD 

-103.208 

44.075 

0.334 

■9 

0.342 

AAM 

-83.656 

42.300 

0.210 

0.148 

0.358 

SHA 

-88.143 

30.694 

0.396 

-0.023 

0.373 

TOL 

-4.049 

39.881 

0.120 

0.300 

0.420 

BOZ 

-111.633 

45.600 

0.046 

0.442 

0.488 

ESK 

-3.205 

55.317 

0.084 

0.415 

0.499 

AQU 

13.403 

42.354 

-0.102 

0.602 

0.500 

1ST 

28.996 

41.046 

0.184 

0.374 

0.558 

MSH 

59.588 

36.311 

0.384 

0.185 

0.569 

HLW 

31.342 

29.858 

-0.256 

0.580 

COP 

12.433 

55.683 

0.174 

0.544 

0.718 
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Figure  5  shows  the  scatter  plot  of  3  different  types  of  station  mb  s  for  Novaya  Zemlya 
explosion  791018.  The  39  good  recordings,  7  noise,  and  14  clips  are  shown  with  filled  circles, 
V-shaped  downward  arrows,  and  upward  arrows,  respectively.  The  raw  station  mb s  (top)  have 
a  standard  deviation  of  0.29  m.u.  Applying  the  “primary”  station  corrections  (cf.  the  “Rev” 
column  in  Table  4;  “S"  term  in  [3])  reduces  the  scatter  to  0.21  m.u.  Applying  the  “secondary" 
corrections  (cf.  the  "Path"  column  in  Table  4;  “F"  term  in  [3])  to  count  for  the  propagation 
effects  reduces  the  scatter  further  down  to  0.1 1  m.u.  The  dashed  lines  of  1  a  range  around 
the  network-averaged  mb  clearly  illustrate  the  remarkable  reduction  of  fluctuation  across  the 
recording  stations.  The  mean  event  mb  itself  is  not  significantly  changed,  however.  Among 
28  Novaya  Zemlya  events  used  in  this  study,  this  event  has  the  smallest  scatter  in  the  result- 
ing  m2.9  values.  The  dramatic  reduction  of  variation  from  m,  to  shows  a  factor  of  nearly 
2.7,  as  compared  to  the  worst  case  of  about  1.26  for  the  event  801011  (Figure  9).  Novaya 
Zemlya  events  typically  exhibit  a  reduction  factor  around  2  (e.g. ,  Figure  7).  Figure  6  shows 
the  scatter  plot  of  3  different  types  of  station  mb  s  for  Novaya  Zemlya  explosion  780927. 
Among  28  Novaya  Zemlya  events  used  in  this  study,  this  event  shows  the  most  dramatic 
reduction  of  variation  with  a  factor  of  2.8. 

Figures  8  and  9  are  the  same  as  Figure  5  except  for  the  events  661027  and  801011, 
respectively.  Although  these  two  events  do  not  show  as  dramatic  reduction  in  variation  as  do 
other  events,  a  factor  of  1 .4  still  illustrates  the  robustness  of  our  joint  inversion  scheme.  The 
seismic  disturbance  of  801011  is  known  to  be  caused  by  two  possibly  simultaneous  explo¬ 
sions  of  7  km  apart  (Lilwall  and  Marshall,  1986;  Stewart  and  Marshall,  1988).  The  resulting 
P- wave  amplitude  pattern  is  significantly  different  from  the  typical  single  shots  from  the  same 
source  region,  which  is  why  the  station/path-correction  procedure  does  not  yield  a  perfor¬ 
mance  as  good  as  in  other  cases.  Note  that  the  path  correction  proposed  in  this  study  not 
.  only  reduces  the  mb  scatter  at  stations  that  reported  the  good  signals,  but  it  also  improves  the 
data  consistency  of  the  censored  recordings,  as  indicated  by  the  shifting  of  the  clipped  record¬ 
ings  (e.g. ,  Figures  5  and  6). 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  791018Z 


Figure  5.  Scatter  plot  of  3  ditferent  types  ot  station  mb s  for  Novaya  Zemlya  explosion  791018.  The  39 
good  recordings,  7  noise,  and  14  clips  are  shown  with  filled  circles,  Y-shaped  downward  arrows,  and 
upward  arrows,  respectively.  The  raw  station  mb s  (top)  have  a  standard  deviation  of  0.29  m.u.  Applying 
the  “primary"  station  corrections  reduces  the  scatter  to  0.21  m.u.  Applying  the  proposed  "secondary" 
corrections  to  count  for  the  path  effects  reduces  the  scatter  further  down  to  0.1 1  m.u.  The  dashed  lines 
around  the  network-averaged  mb  clearly  illustrate  the  remarkable  reduction  of  fluctuation  across  the 
recording  stations.  The  mean  event  mb  itself  is  not  significantly  changed,  however.  Among  28  Novaya 
Zemlya  events  used  in  this  study,  this  event  has  the  smallest  scatter  in  the  resulting  rr^  9  values.  The 
dramatic  reduction  of  variation  from  m,  to  rr^ 9  shows  a  (actor  of  nearly  2  7,  as  compared  to  the  worst 
case  of  about  1.26  for  the  event  801011  (Figure  9).  Novaya  Zemlya  events  typically  exhibit  a  reduction 
factor  around  2. 


VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  780927 Z 
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Figure  6.  Scatter  plot  of  3  different  types  of  station  mb s  for  Novaya  Zemlya  explosion  780927.  Note 
that  the  raw  station  mb s  (top)  have  a  standard  deviation  of  0.32  m.u.,  whereas  the  rr^  g  have  a  standard 
deviation  of  0.11  m.u.  Among  28  Novaya  Zemlya  events  used  in  this  study,  this  event  shows  the  most 
dramatic  reduction  of  variation  with  a  factor  of  2.8. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  87 08 02 Z 


0  SO  100  ISO  200  250  300  350 


Figure  7.  Same  as  Figure  5  except  (or  the  event  870802.  The  variation  reduction  with  a  (actor  o(  2.2  is 
typical  (or  Novaya  Zemlya  explosions. 


-  24- 


VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  661 027 Z 


Raw  Station  rn.'s  (No  Correction  Applied ) 


Azimuth  ( degrees ) 


Figure  8.  Same  as  Figure  5  except  for  the  event  661027.  Although  this  event  does  not  show  as 
dramatic  reduction  in  variation  as  do  other  events,  a  factor  of  1.4  still  illustrates  the  robustness  of  our 
joint  inversion  scheme.  Note  that  the  path  correction  proposed  in  this  study  not  only  reduces  the  mb 
scatter  at  stations  that  reported  the  good  signals,  but  it  also  improves  the  data  consistency  of  the  cen¬ 
sored  recordings,  as  indicated  by  the  shifting  of  the  clipped  recordings. 
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VARIOUS  WWSSN  MAGNITUDES  OF  EVENT  801011Z 


0  50  100  ISO  200  250  300  350 


0  50  100  150  200  250  300  350 

Azimuth  (degrees) 


Figure  9.  Same  as  Figure  5  except  for  the  event  801011.  This  seismic  disturbance  is  Known  to  be 
caused  by  two  possibly  simultaneous  explosions  of  7  km  apart  (Lilwatl  and  Marshall,  1986;  Stewart  and 
Marshall,  1988).  The  resulting  P- wave  amplitude  pattern  is  significantly  different  from  the  typical  single 
shots  from  the  same  source  region,  which  explains  why  the  station/path-correction  procedure  does  not 
yield  a  performance  as  good  as  in  other  cases. 
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5.  YIELD  ESTIMATES  OF  NOVAYA  ZEMLYA  EXPLOSIONS 


Without  calibration  shots,  the  absolute  yields  of  Novaya  Zemlya  explosions  would  have 
to  be  estimated  with  some  assumptions.  Dahlman  and  Israelson  (1977)  used  the  announced 
yields  of  Soviet  PNEs  as  the  baseline.  Burger  etal.  (1986)  estimated  the  yields  by  assuming 
an  average  t’  and  log('P„)  :yield  scaling  relation.  The  yield  estimates  from  Sykes  and  Ruggi 
(1989)  assume  a  mb  bias  of  0.351  m.u.  relative  to  NTS.  Nuttli  (1988)  estimated  the  yields 
using  his  quadric  equation  (Nuttli,  1986a),  which  is  based  on  the  assumption  that  the  Lg: yield 
relationship  is  transportable  to  NNZ  without  any  correction. 

In  this  section,  we  assume  that  the  Pa  scaling  at  Novaya  Zemlya  is  about  the  same  as 
that  at  Semipalatinsk.  For  events  after  1 976,  the  mb  :yield  relationship  at  Matochkin  Shar  of 
Novaya  Zemlya  is  solved  with  that  at  Semipalatinsk  along  with  the  yield-dependent  Pmax/Pa 
ratio  at  these  two  sites.  For  larger  events  (mainly  detonated  before  TTBT  became  effective  in 
1976),  we  combine  the  theoretical  log('F-)  .yield  scaling  in  Bache  (1982)  and  our  mb  : log('F„) 
relationship  using  the  relative  source  strengths  determined  by  Burger  etal.  (1976). 

A  number  of  early  studies  appear  to  support  the  assumption  of  comparable  coupling  and 
mantle  condition.  Geologically,  the  rocks  in  the  upper  crust  of  Novaya  Zemlya  are  of  sedimen¬ 
tary  origin,  metamorphosed  to  lower  greenschist  facies,  that  may  have  physical  properties 
similar  to  those  of  Soviet's  Balapan  test  site  (Leith  etal.,  1990).  DARPA  (1981,  page  32) 
asserts  that  the  mantle  below  Novaya  Zemlya  has  significantly  less  attenuation  as  compared 
to  the  Basin  and  Range  structure.  Analysis  of  modem  digital  array  data  also  favors  a  low  t* 
source  region  (Der  et  ai,  1985). 

Based  on  the  published  yields  (Bocharov  etal.,  1989),  the  mb(Pa): yield  relationship  at 
KTS  is: 

mb(Pa)  *  0.802(±0.020)  log(W)  +  3.834(±0.032)  [6] 

Regressing  m^P^x)  and  mb(Pb)  on  mb(Pg)  with  NNZ  events  after  1976,  we  have:1 

NNZ:  mb(Pb)  =  1 .055(±0.033)  m^Pa)  -  0.056(±0.171)  ,  a  =  0.045  ,  p  =  0.9926  .  [7a] 

NNZ:  m*(Pmax)  =  0.991  (±0.033)  mb{PJ  +  0.446(±0.222)  ,  a  =  0.059  ,  p  =  0.9861  .  [7b] 
Under  the  assumption  that  NNZ  and  KTS  have  similar  coupling  and  mantle  condition  ( i.e . ,  the 

'For  KTS  events  after  1976,  =  0.946(±0.025)  m^P,)  +  0.564(±0.131)  (o  =  0.050,  p  =  0.9908),  and  mdpm*,)  * 

0.923(10.033)  m^Pa)  +  0.912(10.171)  (o  =  0.065,  p  =  0.9839). 
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mb(Pa)  scaling  is  transportable  without  any  correction),  the  mb{P £  :yield  relationship  at  NNZ 
can  be  determined  with  [1]  and  [2a]  as: 

NNZ:  mb(PJ  =  0.846  log(W)  +  3.989  [8a] 

Similarly,  the  mtJiPmax)  yield  relationship  at  NNZ  can  be  determined  with  [6]  and  [7b]  as 

NNZ:  mb(Pmax)  =  0.795  log(W)  +  4.245  [8b] 


Table  5.  Expected  mZ9  for  Various  Test  Sites 

Phase/site 

1  KT 

10  KT 

50  KT 

100  KT 

mb(Pmw)  (NNZ) 

4.245 

5.040 

5.596 

5.835 

mb(Pmax)  (KTS) 

4.399 

5.167 

5.704 

5.935 

mb(P nax)(NTS) 

3.954 

4.771 

5.342 

5.588 

mb(Pb)  (NNZ) 

3.989 

4.835 

5.426 

5.681 

mb(Pb)  (KTS) 

4.130 

4.930 

5.489 

5.730 

mb(Pb)(  NTS) 

3.674 

4.505 

5.086 

5.336 

mb(Pa)  (KTS, NNZ) 

3.834 

4.636 

5.197 

5.438 

mb(Pa)(  NTS) 

3.607 

4.372 

4.907 

5.137 

Table  5  lists  the  expected  mb  values  for  each  of  Pa,  Pb,  and  Pmax  phases  from  NTS, 
NNZ,  and  KTS  explosions  based  on  Equations  [6]  and  [8].  The  estimated  “mean”  m2  9  bias 
can  then  be  computed  in  a  straightforward  manner  (Table  6).  The  bias  estimates  based  on 
m22  can  be  found  in  Jih  and  Wagner  (1991b).  Leith  etal.  (1990)  suggest  that  the  NNZ 
explosions  are  primarily  emplaced  in  near-horizontal  tunnels,  excavated  into  the  steep  moun¬ 
tain  slopes.  Thus  it  would  seem  plausible  that  the  difference  in  the  free-surface  interaction  at 
NNZ  and  Balapan/Murzhik  could  contribute  much  of  the  mb  bias  in  Pmax  and  Pb  we  are 
observing.  The  majority  of  NNZ  explosions  after  1976  are  clustered  at  yields  between  50  and 
70  KT  (cf.  Table  8).  Within  this  range  of  yield,  the  expected  mb  bias  is  about  -0.1 1  m.u.  rela¬ 
tive  to  KTS  using  mtj(iPm „*),  which  can  explain  satisfactorily  why  Evernden  and  Marsh  (1987) 
and  Sykes  and  Cifuentes  (1984)  observe  the  same  mb  bias  of  about  0.1  m.u.  between  Semi- 
palatinsk  and  Novaya  Zemlya  on  their  Ms:mb  plot. 
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Table  6.  Mean  /r^  9  Bias 

Phase 

1  KT 

10  KT 

50  KT 

100  KT 

mb(Pmax)  (KTS-NTS) 

0.45 

0.40 

0.36 

0.35 

mb(Pb)(KTS-HTS) 

0.46 

0.43 

0.40 

0.39 

mb(PJ{  KTS-NTS) 

0.23 

0.26 

0.29 

0.30 

m6(Pmax)(  NNZ-NTS) 

0.29 

0.27 

0.25 

0.25 

m^P*,)  (NNZ-NTS) 

0.32 

0.33 

0.34 

0.35 

mb(Pa)  (NNZ-NTS) 

0.23 

0.26 

0.29 

0.30 

Burger  et  al.  (1 986)  estimated  the  relative  source  strengths  of  1 1  NNZ  events  detonated 
before  1976  based  on  the  intercorrelation  method  with  the  event  671021  as  a  reference 
(Table  7).  Regressing  our  path-corrected  mb  values  on  these  relative  source  strengths,  we 
have: 

mb(Pa)  =  0.854{±0.048)  log['F../4'„  (671 021 )]  +  5.472(±0.032)  [9a] 

mb[Pb)  =  0.900(±0.048)  log[^„/'F.(671021)]  +  5.672(±0.033)  [9!  ] 

mb(Pmax)  =  0.856(10.047)  logJ'P. /¥.  (671021)]  +  5.856(10.032)  [9c] 
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|  Table  7.  Relative  Source  Size  'P-/H,oo(671021)  of  Novaya  Zemlya  Explosions*  j 

Event 

'f'«./'t*»(671021) 

S.E.M. 

661027 

5.26 

0.43 

671021 

1.00 

0.11 

681107 

1.53 

0.12 

691014 

2.22 

0.16 

701014 

11.29 

1.01 

710927 

7.47 

0.66 

720828 

4.10 

0.36 

730912 

16.24 

2.20 

740829 

5.44 

0.48 

750823 

5.28 

0.45 

751021 

4.96 

— 

0.56 

*:  from  Burger  at  al. .  (1986). 


In  deriving  Equation  [9],  the  standard  errors  in  both  mb  and  VF„  measurements  (shown 
in  Tables  1  and  7,  respectively)  are  taken  into  account  using  a  generalized  regression  pro¬ 
cedure  described  in  Jih  and  Wagner  (1991b)  and  Jih  etal.  (1991).  Combining  the  theoretical 
logt'PJ  .yield  scaling  of  log('PJ  =  0.89  log(W)  +  C  (e.g. ,  Bache,  1982),  we  can  infer  the 
mb  .yield  relationship  for  larger  Novaya  Zemlya  explosions  as 

mb(Pa )  =  0.760  log(W)  +  3.920  [10a] 

mb(Pb)  =  0.801  log(W)  +  4.079  [10b] 

mb(Pmax)  =  0.762  log(W)  +  4.312  [10c] 

Equation  [10a]  is  used  if  mb(Pa)  is  larger  than  5.472.  Similar  turning  points  for  mb(Pb)  and 
mb{Pmax)  are  set  at  5.672  and  5.856,  respectively,  which  are  the  corresponding  intercepts  in 
Equations  [9a]  through  [9c].  At  these  turning  points,  Equations  [8]  and  [10]  give  identical 

result,  as  we  have  used  these  turning  points  to  constraint  the  unknown  constant  C  in  the 

theoretical  logf'PJ  scaling. 

Figures  10  through  12  show  the  regression  of  /TT29(Pa),  m29(Pb),  and  m:n(Pma%)  on  the 
relative  source  size  4'„/vF<»(671021)  determined  by  Burger  et  al.  (1986),  which  correspond  to 
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Equations  [9a]  through  [9c],  respectively.  The  uncertainties  in  the  mb  s  and  the 
vf'«/'Fo.(671021)  are  taken  into  account  through  800  bootstrap  resamplings.  The  darkened 
bundle  is  actually  the  collection  of  all  800  regressions,  each  produced  by  a  possible  realiza¬ 
tion  of  11  perturbed  ( mb ,  (671021))  pairs.  The  95%  confidence  band  (shown  as  2 

curves  around  the  darkened  bundle)  is  narrower  near  the  centroid  and  wider  towards  both 
ends,  as  expected.  The  individual  95%  confidence  intervals  of  the  two  inferred  parameters 
(i.e. ,  the  slope  and  the  intercept  of  the  calibration  curve)  are  shown  with  the  dashed  line  in 
the  scatter  plot  (bottom).  Note  that  the  dashed  rectangle  is  not  the  joint  90%  confidence  inter¬ 
val,  however,  due  to  the  highly  correlated  nature  of  the  two  parameters.  The  Pb  phase 
appears  to  fit  the  relative  source  sizes  slightly  better  than  do  Pa  and  ^rtiax  • 

Table  8  gives  our  yield  estimates  based  on  the  m2  9  measurements  listed  in  Table  1.  For 
comparison,  the  yields  estimated  by  Nuttli  (1988),  Sykes  and  Ruggi  (1989),  and  Burger  et  al. 
(1986)  are  also  included.  Most  of  Nuttli’s  (1988)  yield  estimates  based  on  his  quadric  formula 
tend  to  be  larger,  as  noted  by  Lay  (1991).  Sykes  and  his  colleagues  (  Sykes  and  Wiggins, 
1986;  Sykes  and  Davis,  1987;  Sykes  and  Ruggi,  1989)  examined  the  apparent  clustering  of 
certain  yields  of  Novaya  Zemlya  explosions  in  the  context  of  particular  warhead  types,  and 
they  estimated  the  peak  value  as  500  and  90  KT  for  events  before  and  after  1976,  respec¬ 
tively.  Our  peak  value  for  recent  events  lies  between  Sykes’  90  KT  and  Israelson’s  (1991)  50 
KT.  To  lower  down  our  estimate  of  the  peak  value  would  imply  that  NNZ  has  basic  coupling 
much  more  efficient  than  that  at  KTS  or  that  NNZ  has  a  less  attenuating  mantle  than  KTS, 
which  would  seem  unlikely.  Therefore  Israelson’s  (1991)  peak  value  of  50  KT  could  be  some¬ 
what  low.  Recent  laboratory  experiments  by  Miller  and  Florence  (1991)  indicate  that 
saturated  limestones  under  the  frozen  condition  could  have  a  poorer  coupling.  Thus  there  is  a 
possibility  that  yield  estimates  of  Novaya  Zemlya  events  based  on  every  seismic  method 
could  be  slightly  underestimated,  and  all  the  peak  values  might  need  be  increased  somewhat 
(Blandford,  personal  communication). 

There  is  a  clustering  of  four  NNZ  events  around  4,„/'P„(671021)  =  5.235  before  1976 
(Figures  10  through  12),  which  corresponds  to  577,  510,  and  556  KT  in  Equations  [10a], 
[10b],  and  [10c],  respectively.  Between  1973  and  1976  the  U.S.S.R.  deployed  five  strategic 
systems  with  warhead  yields  between  300  and  600KT  according  to  some  non-seismic  source 
of  information  (Samuel,  1985;  Sykes,  1985;  Evernden  and  Marsh,  1987),  which  provides  an 
indirect  support  of  our  estimate  of  550KT  and  Sykes’  500KT  for  the  repeated  tests  at  NNZ 
prior  to  TTBT  era. 
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mb:Yield  Relationship  at  Novaya  Zemlya 


DWLS  (uncertain  X  &  Y):  s=0.85(0.048),  1=5.47(0.032),  11.  data  used, 
95%  error  tn  Y  at  X=1 ,5,10,50,100:  0.18, 0.12,  0.15,  042, 0.40, 

95%  factor  In  X  at  X=1 ,5,10,50,1 00:  2.69, 1 .94, 249, 5.67, 8.70 

OWLS  (precise  X  assumed):  S=047(0.051),  1=5.46(0436) 

Standard  LS:  S  =0.86(0.056),  1=5.47(0.041) 


Scatter  Plot  of  Inferred  Parameters 


95%  confidence  Interval  of  slope:  0.85+/-0.10B 
95%  confidence  Interval  ot  Intercept:  5.47+/-0.073 
[97.5%  quantile  of  t(9.  D.o.F.),  2.262,  used] 


Figure  10.  Regressing  the  flk9(Pa)  °n  the  relative  source  size  4,>/4/„(67l02l)  determined  by  Burger 
et  al.  (1986).  The  uncertainties  in  the  mb s  and  the  VF>/VP-(671021)  are  taken  into  account  through  800 
bootstrap  resamplings.  The  darkened  bundle  is  actually  the  collection  of  all  800  regressions,  each  pro¬ 
duced  by  a  possible  realization  of  11  perturbed  (mb,  VF>/'P_( 671021))  pairs.  The  95%  confidence  band 
(shown  as  2  curves  around  the  darkened  bundle)  is  narrower  near  the  centroid  and  wider  towards  both 
ends,  as  expected.  The  individual  95%  confidence  intervals  of  the  two  inferred  parameters  ( i.e . ,  the 
slope  and  the  intercept  of  the  calibration  curve)  are  shown  with  the  dashed  line  in  the  scatter  plot  (bot¬ 
tom).  Note  that  the  dashed  rectangle  is  not  the  joint  90%  confidence  interval,  however,  due  to  the 
highly  correlated  nature  of  the  two  parameters. 
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mb:Yield  Relationship  at  Novaya  Zemlya 


OWLS  (uncertain  X  &  Y):  S=0.90(0.048),  1=5.67(0.033),  1 1 .  data  used, 
95%  error  In  Y  at  X=1 ,5,10,50,100:  0.17, 0.12, 0.14, 0.29, 0.36, 

95%  factor  In  X  at  X=1 ,5,10,50,100:  2.37, 1 .83, 2.0S,  4.46, 6.37 

OWLS  (precise  X  assumed):  S=0.91  (0.038),  1=5.66(0.027) 

Standard  LS:  S=0.91  (0.042),  1=5.67(0.031) 


Scatter  Plot  of  Inferred  Parameters 


95%  confidence  Interval  of  slope:  0.90+/-0.108 
95%  confidence  Interval  of  Intercept:  5.67+/-0.074 
(97.5%  quantile  of  t{9.  D.o.F.),  2.262,  used] 


Figure  11.  Same  as  Figure  10  except  that  /7Vg(Pt,)  are  used.  The  Pb  phase  appears  to  fit  the  relative 
source  sizes  slightly  better  than  do  Pa  and  ^max  • 
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800-bootstrap  Intercepts  mb(WWSSN,ML8,Pmax) 


mb:Yietd  Relationship  at  Novaya  Zemlya 


OWLS  (uncertain  X  &  Y):  8=0.86(0.047),  1=5-86(0.032),  11 .  data  used, 
95%  error  In  Y  at  X=1 ,5,1 0,50,100:  0.19, 0.12, 0.16, 0.33,  0.41, 

95%  factor  in  X  at  X=1 ,5,10,50,100:  2.74,1.94, 2.31, 5.75, 8.85 

OWLS  (precise  X  assumed):  S=0.87(0.053),  1=5.84(0.038) 

Standard  LS:  S=0.86(0.058),  1=5.85(0.043) 


Scatter  Plot  of  Inferred  Parameters 
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800-bootstrap  Slopes 


95%  confidence  Interval  of  slope:  0.86+/-0.106 
95%  confidence  Interval  of  Intercept:  5.86+/-0.072 
[97.5%  quantile  of  t(9.  D.o.F.),  2.262,  used] 


Figure  12.  Same  as  Figure  10  except  the  TThs(Pmiu)  are  used. 
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Event 


Table  8.  Yield  Estimates  of 


Yield  [KT]  (this  study) 


Novaya 


Date 

Pa 

Pb 

P 

^max 

661027 

673 

602 

642 

671021 

96 

82 

86 

681107 

164 

161 

180 

691014 

270 

225 

245 

701014 

2058 

1537 

1888 

710927 

1140 

893 

945 

720828 

518 

488 

498 

730912 

1639 

1764 

1693 

740829 

834 

760 

917 

750823 

862 

766 

760 

751021 

748 

647 

820 

760929 

58 

56 

51 

761020 

3 

4 

5 

770901 

46 

55 

50 

771009 

3 

2 

2 

780810 

93 

89 

109 

780927 

38 

42 

41 

790924 

66 

71 

76 

791018 

68 

61 

67 

801011 

49 

52 

61 

811001 

63 

63 

63 

821011 

40 

35 

32 

830818 

70 

65 

69 

830925 

58 

56 

57 

841025 

49 

55 

54 

870802 

71 

66 

65 

880507 

42 

37 

34 

881204 

63 

70 

63 

Zemlya  Explosions 


Yield  [KT]  (earlier  study) 


SR1  mb(La)2  BBL3 


1:  from  Sykes  and  Ruggi  (1989);  based  on  body-wave  magnitudes. 
2:  from  Nuttli  (1988);  based  on  the  quadric  fit. 

3;  from  Burger  et  at.  (1986);  based  on  Amchitka  scaling. 
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5.  YIELD  ESTIMATES  WITH  NUTTLI’S  L„  MEASUREMENTS 

fl 

Lay  (1991)  suggests  that  there  still  exists  significant  uncertainty  in  yield  estimates  for 
many  of  the  early  large  events  at  Novaya  Zemlya.  We  note  that  Nuttli’s  (1 988)  yield  estimates 
are  systematically  larger  than  those  inferred  by  Burger  etal.  (1986),  while  our  yield  estimates 
for  the  larger  historical  events  seem  to  lie  between  these  two  (with  730912  as  an  exception; 
cf.  Table  10).  A  few  experiments  with  Nuttli’s  (1988)  m^Lg)  measurements  have  been  con¬ 
ducted  to  explore  some  of  the  possibilities. 

The  first  question  to  be  examined  is  whether  Nuttli’s  (1986a)  Lg  scaling  formulae  are 
adequate.  Jih  etal.  (1990)  point  out  that  it  is  not  quite  clear  how  Nuttli's  (1986a)  derived  his 
original  linear  calibration  formula.  Nuttli's  data  set  of  high-coupling  NTS  shots  (cf.  Nuttli, 
1986a,  page  2144)  included  the  Pahute  Mesa  event  HANDLEY  which  has  an  announced 
yield  of  >1000KT.  However,  Nuttli  seemed  to  have  treated  the  yield  as  exactly  1000KT  in  his 
calculations  (cf.  Figures  7  and  9  of  Nuttli,  1986a).  Jih  etal.  (1990)  tested  eight  possible  com¬ 
binations  with  Nuttli’s  mb(Lg)  measurements: 

•  including  the  granite  events  (SHOAL  and  PILEDRIVER)  or  not, 

•  limiting  m£Lg)  to  [5.2, 6.7]  or  not, 

•  assuming  HANDLEY  was  1000KT  or  deleting  HANDLEY  from  the  regression. 

None  of  the  eight  experiments  could  give  an  exactly  identical  formula  to  the  linear  fit  given  by 
Nuttli  (1986a).  It  seems  very  likely  that  Nuttli  might  have  regressed  the  published  yields  on 
his  mb(Lg)  values  (i.e.,  the  so-called  “ /-regression”  models)  with  some  unspecified  constraint 
on  the  data  set. 

A  more  extensive  data  set  of  NTS  mb(Lg)  values  has  been  compiled  by  R.  Geil  at 
Lawrence  Livermore  National  Laboratory  [LLNL]  (Patton,  1988).  This  data  set  consists  of  47 
events  below  the  water  table  on  Pahute  Mesa  and  Yucca  Flat  of  NTS.  Patton  (1988) 
regressed  this  data  set  against  the  official  yields,  and  the  calibration  curve  meets  that  of 
Nuttli's  (1986a)  at  mb(Lg)  =6.028  or  178KT.  Beyond  that  level,  yield  estimates  based  on  Geil’s 
data  set  are  even  larger  than  that  by  Nuttli  (1986a).  Patton  (1988)  repeated  Nuttli’s  (1986a) 
procedure  to  estimate  the  yields  of  69  NTS  high-coupling  shots  using  RMS  Lg  magnitudes 
measured  at  LLNL’s  high-quality  regional  digital  network.  The  resulting  scaling  curve  is  nearly 
parallel  to  that  of  Nuttli’s,  and  hence  yield  estimates  based  on  Patton’s  formula  are  systemati¬ 
cally  smaller  than  those  with  Nuttli’s  or  Geil's  results  (Tables  9  and  10).  For  instance,  based 
on  Patton’s  regression  result,  the  predicted  mb(Lg)  at  explosive  yields  of  10,  50,  100,  150KT 
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are  5.159,  5.687,  5.914,  and  6.047,  respectively.  Nuttli’s  (1986a)  original  regression  with  22 
NTS  shots  recorded  at  WWSSN  stations  gives  5.072,  5.607,  5.837,  and  5.972,  respectively. 
The  NTS  explosions  Patton  used  are  clustered  around  mb(Lg)  =  5.8.  Beyond  that  level,  the 
difference  in  yield  estimates  between  Nuttli's  and  Patton’s  predictions  are  by  no  means  negli¬ 
gible.  The  data  recorded  at  LLNL’s  regional  digital  network  probably  have  a  quality  better 
than  those  WWSSN  film  chips  which  Nuttli  (1986a)  read.  Perhaps  it  is  reasonable  to  expect 
that  the  Patton’s  (1988)  mb(Lg)  :yield  calibration  curve  is  more  accurate  than  Nuttli’s.  Table  9 
gives  the  predicted  yield  at  several  mb(Lg)  levels,  with  the  applicable  range  of  each  scaling 
formula  ignored. 


Table  9.  Expected  Yields  [KT]  at  Given  mb(Lg)  Values 

Scaling 

4.5 

5.0 

5.5 

o 

<6 

6.5 

NuttiKigaea)1 

3 

10 

37 

152 

779 

Nuttli(1986a)2 

2* 

8* 

36 

163 

736 

Geil3 

1* 

5* 

29 

162 

885 

Patton(1988)4 

1 

6 

28 

130 

597 

1)  Nuttli’s  (1986a)  quadric  formula:  m£La)  •  3.943  +  1.124  log(W)  -  0.0829  [tag(W)f. 

2)  Nuttli  (1986a):  m^Lg)  -  4.307  [±0.067]  +  0.765  [±0.027]log(W)  for  5.2  <  m^Lg)  <  6.7. 

3)  Gail  (Patton.  1988):  m^Lg)  -  4.505  [±0.067]  +  0.677  [±0.029]k>g(W)  for  5.28  <  m^Lg)  <  6.65. 

4)  Patton  (1988):  m^Lg)  -  4.404  [±0.048]  +  0.755  [±0.022]log(W)  for  4.22  <  m^Lg)  <  6.7. 

*)  Extrapolated  value. 


Table  10  gives  the  yield  estimate  of  NNZ  events  inferred  from  Nuttli’s  (1988)  m^Lg) 
values  and  various  Lg  scaling  relationship,  assuming  that  these  Lg  scaling  formulae  are  port¬ 
able  and  that  Nuttli’s  mb(Lg)  measurements  of  NNZ  events  are  accurate.  In  comparing  Table 
8  against  10,  we  do  see  some  indication  that  substituting  Nuttli’s  (1988)  NNZ  m^Lg)  values 
into  Patton’s  (1988)  Lg  scaling  for  NTS  shots  could  give  yield  estimates  in  better  agreement 
with  our  yields,  particularly  with  those  based  on  the  Pb  phase.  However,  further  investigation 
is  definitely  necessary. 
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Table  10.  Yield  Estimates  Based  on  Nuttli’s  m^Lg)  Measurements 


Event 

Yield  [KT] 

Date 

m^LQ) 

Nuttli1 

Nuttli2 

Nuttli3 

Geil4 

Patton5 

640918 

4.37 

2.5 

2 

1* 

1* 

1 

641025 

5.19 

16.4 

17 

14 

10* 

11 

661027 

6.45 

644 

653 

633 

746 

513 

671021 

6.06 

180 

182 

196 

198 

156 

681107 

6.17 

253 

257 

272 

288 

218 

691014 

6.31 

399 

405 

415 

464 

335 

701014 

6.75 

1970 

2000 

1561 

2071 

1280 

710927 

6.68 

1500 

1519 

1265 

1632 

1034 

720828 

6.42 

580 

588 

578 

674 

468 

730912 

6.89 

3510 

3572 

2379* 

3333* 

1962* 

740829 

6.60 

1110 

1122 

994 

1243 

810 

750823 

6.47 

690 

700 

672 

799 

545 

751021 

6.43 

600 

609 

596 

697 

482 

760929 

5.83 

91 

92 

98 

91 

77 

761020 

5.24 

19 

19 

17 

12* 

13 

770901 

5.93 

122 

123 

132 

127 

105 

771009 

4.99 

10 

10 

8 

5 

6 

780810 

5.83 

91 

92 

98 

91 

77 

780927 

5.69 

61 

62 

64 

56 

51 

790924 

5.79 

81 

82 

87 

79 

69 

791018 

5.78 

79 

80 

84 

76 

66 

801011 

5.77 

76 

77 

82 

74 

64 

811001 

5.91 

116 

116 

125 

119 

99 

821011 

5.78 

79 

80 

84 

76 

66 

830818 

5.99 

145 

147 

158 

156 

126 

830925 

5.86 

99 

100 

107 

100 

85 

1)  Nuttli-furnished  yields  based  on  his  quadric  formula  (page  881  of  Nuttli,  1988). 

2)  Recomputed  estimates  with  Nuttli's  (1986a)  quadric  formula  (see  footnote  of  Table  9). 

3)  Nuttli's  (1986a)  linear  fit  (see  footnote  of  Table  9). 

4)  Geil  (Patton,  1988)  (see  footnote  of  Table  9). 

5)  Patton  (1988)  (see  footnote  of  Table  9). 

*)  Extrapolated  value. 
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Lay  (1991)  compared  a  megaton-level  Amchitka  event  MILROW  with  two  Novaya  Zemlya 
events  of  comparable  Ms ,2  and  asserts  that  there  is  significant  discrepancy  in  the  yield  esti¬ 
mates  based  on  body  wave  and  surface  wave.  We  feel  a  little  more  optimistic  on  this  issue,  if 
the  uncertainty  in  each  estimate  is  taken  into  account.  Table  11  lists  our  m2. 9  of  these  two 
NNZ  events  (from  Table  1)  and  MILROW  (Jih  and  Wagner,  1992).  The  yield  estimates  of 
Burger  et  at.  (1986)  are  in  reasonable  agreement  with  our  estimate  based  on  Pa  (or  Pb, 
Pmax’  cf-  Table  8).  Nuttli’s  (1988)  mb(Lg)  measurements  simply  suggest  these  two  NNZ 
events  could  have  very  comparable  size  in  terms  of  Lg  excitation  alone. 


Table  1 1 .  Comparison  of  MILROW  and  2  Novaya  Zemlya  Explosions 

Event 

n^g  (this  study) 

Yield  Estimate 

(Pmax  ) 

mb(Pb) 

mb(Pa) 

Pa 

BBL 

Geil 

Patton 

MILROW 

6.618 

6.321 

6.071 

_ 

_ 

701014 

6.808 

6.631 

6.438 

2058 

1714 

2071 

1280 

710927 

6.579 

6.442 

6.243 

1140 

973 

1632 

1034 

zMarshaJt  et  al.  (1979)  and  Douglas  et  a!.  (1987)  give  Ms  of  MILROW  as  5.05  and  5.2,  respectively.  Sykes  and  Ruggi 
(1989)  estimate  the  Ms  of  701014Z  and  710927Z  as  5.02  and  5.06.  respectively. 
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6.  CONCLUSIONS  ANO  RECOMMENDATIONS 


Along  with  an  extensive  data  set  of  worldwide  explosions  recorded  at  a  global  network, 
teleseismic  body-wave  amplitudes  from  28  Novaya  Zemlya  explosions  are  measured  and 
analyzed  to  isolate  the  propagation  characteristics  and  to  derive  a  better  measure  of  the 
source  size.  This  new  mb  factoring  procedure  provides  more  stable  mb  measurements 
across  the  whole  recording  network  with  a  reduction  in  the  fluctuational  variation  by  a  factor  of 
up  to  3.  The  variation  reduction  is  typically  a  factor  of  2,  except  for  a  known  double  explosion 
which  has  a  factor  of  1.26.  In  principle,  this  procedure  can  be  applied  to  other  types  of 
network-recorded  magnitudes  as  well,  such  as  mb  (the  band-passed  spectral  amplitude,  see 
Bache,  1982,  and  Murphy  eta!.,  1989),  m^Lg),  M0,  and  Ms. 

Our  result  indicates  that  paths  from  the  northern  test  site  in  Novaya  Zemlya  to  stations  in 
North  America  have  systematically  faster  arrivals  and  smaller  amplitudes,  suggesting  a  pro¬ 
found  defocusing  effect  on  the  first  arrivals;  while  stations  in  Ireland,  Scotland,  Spain,  Ban¬ 
gladesh,  northern  India,  Pakistan,  Korea,  and  Kenya  report  slow  arrivals  and  large  amplitudes, 
suggesting  a  focusing  effect.  Amplitudes  for  paths  to  Greenland,  Iceland,  Alaska,  Turkey,  Ger¬ 
many,  Luzon,  Zimbabwe,  Italy,  Pueto  Rico,  Ethiopia,  and  Hawaii,  however,  seem  to  be  con¬ 
trolled  by  the  anelastic  attenuation  with  slow  rays  also  associated  with  small  amplitudes,  and 
fast  rays  associated  with  large  amplitudes.  For  paths  showing  strong  attenuation  or  strong 
defocusing,  our  observations  favor  the  “weak  signal”  hypothesis  which  predicts  complex 
waveforms  for  such  paths.  The  separation  of  path  effects  from  station  effects  also  provides 
direct  clues  as  how  an  old  fold-belt  structures  (like  Novaya  Zemlya)  could  "modulate"  the 
short-period  P- wave  amplitude  and  the  travel-time  patterns,  in  a  manner  similar  to  that  tec¬ 
tonic  release  affects  Ms-  A  strong  correlation  between  P- wave  amplitude  and  Lg  detection  at 
teleseismic  distance  is  also  observed.  Thus  a  thorough  assessment  of  WWSSN’s  remote 
monitoring  capability  along  this  line  could  provide  useful  clues  such  as  which  stations  should 
be  searched  for  the  Lg  phase.  Naturally,  such  "detection  map”  will  be  varying  from  one  test 
site  to  another,  as  discussed  in  Jih  and  Wagner  (1992). 

Although  our  results  can  explain  some  of  the  propagation  complexities  in  the  initial 
P- wave  arrivals,  a  follow-up  study  is  needed  to  quantify  further  the  contribution  of  near¬ 
source  scattering  to  the  waveform  complexity  in  the  P  coda  (such  as  that  in  Lay  and  Welc, 
1987)  which  is  not  covered  in  this  study.  Our  previous  modeling  effort  (McLaughlin  and  Jih, 
1986,  1987,  1988;  Jih  and  McLaughlin,  1988)  of  utilizing  the  linear  finite-difference  code  (Jih 
eta!.,  1988)  focused  on  the  effects  of  mountainous  topography  and  hypothetical 
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heterogeneity  in  the  upper  crust  on  teleseismic  and  regional  phases  with  somewhat  simplied 
structures  of  other  test  sites.  We  suggest  that  the  follow-up  research  be  accompanied  with 
some  well-constrained  forward  modeling  study  using  more  realistic  structures  of  Novaya  Zem- 
lya. 

Assuming  the  basic  coupling  and  the  mantle  condition  at  Novaya  Zemlya  are  comparable 
to  those  at  Eastern  Kazakhstan,  the  yield  estimates  for  the  28  Novaya  Zemlya  events  based 
on  the  path-corrected  mb  presented  in  this  study  range  from  2  to  2100  KT,  with  peak  values 
at  550KT  and  65KT  for  events  before  and  after  1 976,  respectively.  The  relative  source  size 
determined  by  Burger  etal.  (1986)  and  the  theoretical  'f'«:yield  scaling  are  combined  to 
extrapolate  the  mb  scaling  to  higher  end,  and  the  resulting  yield  estimates  are  in  reasonable 
agreement  with  those  in  previous  studies. 

There  appears  to  be  a  bias  of  0.11  magnitude  units  in  mb(Pmax)  -mb(PJ  between 
Eastern  Kazakh  and  Novaya  Zemlya.  This  bias  could  be  largely  due  to  the  difference  in  pP 
interference  at  these  two  test  sites.  Previous  studies  suggest  a  mb  bias  of  about  0.35  m.u. 
between  NNZ  and  NTS,  based  on  the  spectral  slope  study  ( e.g .,  Der  etal.,  1985)  and  the 
mb:Ms  shift  (e.g.,  Evemden  and  Marsh,  1987).  Result  presented  in  this  study  gives  a 
mb(Pmax)  bias  of  about  0.25  m.u.  (NNZ-NTS)  and  0.36  m.u.  (KTS-NTS)  at  50  KT  level.  It 
seems  that  Nuttli’s  (1988)  Balapan-NNZ  bias  (relative  to  m^Lg) )  may  be  slightly  biased  high 
by  about  0.05  m.u.,  which  is  very  likely  to  be  resulted  from  the  inherent  inconsistency  associ¬ 
ated  with  the  ISC  bulletin  mb  values  that  Nuttli  (1986b,  1987,  1988)  used.  Note  that  ISC  bul¬ 
letin  mb  is  just  the  simple  network  average  of  the  raw  station  mb  ( i.e .  in  Equation  [1]) 
values  without  any  further  processing.  Our  study  illustrates  the  advantages  and  importance  of 
adopting  some  sophisticated  post-processing  methodology  in  determining  the  optimal  network 
mb  values,  even  when  the  network  (such  as  WWSSN  and  ISC)  already  possess  a  broad  spa¬ 
tial  coverage. 
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APPENDIX:  YIELD  ESTIMATES  OF  SEMIPALATINSK  EXPLOSIONS 

The  accuracy  of  our  absolute  yield  estimates  of  Novaya  Zemlya  explosions  presented  in 
this  study  rely  on  that  of  Semipalatinsk  explosions  as  well  as  the  validity  of  the  comparable¬ 
coupling  assumption.  It  is  fortuitous  to  have  the  source  information  by  Bocharov  etal.  (1989) 
(and  Vergino,  1989)  to  calibrate  the  Semipalatinsk  test  site.  The  small  scatter  around  the  fol¬ 
lowing  calibration  curves  based  on  the  regression  of  our  path-corrected  mb  on  the  published 
yields  illustrates  how  well  the  fit  can  be  at  the  Central  Asian  test  site.  Detailed  discussion  can 
be  found  in  Jih  and  Wagner  (1992). 

mb{Pa)  =  0.802(±0.020)  log(W)  +  3.834(±0.032) 
mb(Pb)  =  0.800(±0.020)  log(W)  +  4.130(±0.031) 
mb(pmax)  =  0.768(±0.019)  log(W)  +  4.399(10.030) 

We  have  utilized  these  calibration  curves  to  estimate  the  yield  of  all  Semipalatinsk  explo¬ 
sions  in  our  data  set,  and  the  result  is  summarized  in  the  following  Tables.  For  the  cratering 
events  (such  as  6501 15B)  the  yield  estimate  based  on  the  first  motion  (i.e.,  Pa)  should  be 
used,  since  no  depth  correction  (e.g. ,  Marshall  etal. ,  1979)  has  been  applied  to  mb(Pb)  or 
mb(p max)  in  our  study.  Note  the  excellent  match  between  the  announced  yield  (100-150  KT) 
of  6501 15B  and  the  predicted  value  (1 1 1  KT).  Another  example  is  the  Soviet  Joint  Verification 
Experiment  [JVE]  of  88091 4B,  which  is  said  to  have  an  on-site  hydrodynamic  yield  measure¬ 
ment  of  115  KT  (Gordan,  1988).  Both  examples  illustrate,  again,  that  Pa  from  hard-rock  test 
sites  in  stable  region  could  be  a  very  favorable  phase  for  the  source  size  determination.  The 
mb(pmax)  of  25  Balapan  events  for  which  the  mb(Pb)  and  mb(Pa)  are  missing  are  based  on 
mb  values  published  by  Lilwall  etal.  (1988).  A  correction  to  convert  the  B(A)  (cf.  Equation 
[1])  of  Gutenberg  and  Richter  (1956)  to  that  of  Veith  and  Clawson  (1972)  is  applied  to  every 
station  mb  before  these  recordings  are  incorporated  into  our  data  set. 
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Magnitudes  of  Semipalatinsk  Explosions 

Event  #  of  Signals  Magnitudes  (Jih  and  Wagner,  1992)  Yield 

Date  Site  Ns  Nn  Nc  S.E.M.  Pa  Pb  Pmax 


6501 15B 

BTZ 

45  1  2 

0.028 

5.473 

5.709 

5.865 

100-150 

651121D 

Deg 

48  15  1 

0.024 

4.962 

5.240 

5.452 

29 

66021 3D 

Deg 

51  4  10 

0.024 

5.717 

5.965 

6.152 

125 

660320D 

Deg 

49  9  8 

0.024 

5.416 

5.697 

5.916 

100 

660507 D 

Deg 

9  26  1 

0.033 

4.089 

4.237 

4.529 

4 

661019D 

Deg 

51  10  5 

0.024 

5.164 

5.423 

5.596 

20-150 

661218M 

Mzk 

55  8  1 

0.024 

5.395 

5.632 

5.852 

20-150 

670226 D 

Deg 

48  9  6 

0.025 

5.438 

5.688 

5.914 

20-150 

67091 6M 

Mzk 

36  29  2 

0.024 

4.657 

4.937 

5.182 

<20 

670922M 

Mzk 

35  31  1 

0.024 

4.516 

4.840 

5.118 

10 

671122M 

Mzk 

7  63  0 

0.023 

3.975 

4.353 

<20 

68061 9B 

BNE 

28  3  2 

0.034 

4.666 

5.002 

5.256 

<20 

680929 D 

Deg 

50  7  6 

0.025 

5.222 

5.511 

5.710 

60 

690531 M 

Mzk 

30  30  0 

0.025 

4.468 

4.885 

5.115 

<20 

690723D 

Deg 

38  20  1 

0.025 

4.711 

5.022 

5.248 

16 

69091  ID 

Deg 

19  38  0 

0.026 

4.141 

4.381 

4.709 

<20 

691130B 

BTZ 

49  0  0 

0.028 

5.362 

5.733 

5.915 

125 

691228M 

Mzk 

45  9  3 

0.026 

5.264 

5.551 

5.753 

46 

700721 M 

Mzk 

38  20  1 

0.025 

4.689 

5.033 

5.281 

<20 

701104M 

Mzk 

38  22  1 

0.025 

4.934 

5.137 

5.349 

<20 

71 0322 D 

Deg 

43  143 

0.025 

5.117 

5.408 

5.587 

20-150 

710425D 

Deg 

37  5  0 

0.030 

5.434 

5.696 

5.891 

90 

710606M 

Mzk 

38  12  2 

0.027 

4.879 

5.218 

5.425 

16 

Magnitudes  of  Semipalatinsk  Explosions 


Event 

#  of  Signals 

Magniti 

Date 

Site 

Ns  Nn  Nc 

S.E.M. 

710619M 

Mzk 

41  13  0 

0.027 

710630B 

BTZ 

31  19  1 

0.027 

711009M 

Mzk 

27  12  3 

0.030 

71 1021 M 

Mzk 

32  9  0 

0.031 

711230D 

Deg 

1630 

0.045 

72021 OB 

BNE 

34  8  2 

0.029 

720328D 

Deg 

28  17  0 

0.029 

720816D 

Deg 

23  23  1 

0.029 

720826 M 

Mzk 

29  15  2 

0.029 

720902 M 

Mzk 

15  29  0 

0.029 

721102B 

BSW 

42  1  15 

0.026 

721210D 

Deg 

30  7  5 

0.030 

721210B 

BNE 

44  2  11 

0.026 

730723B 

BTZ 

52  1  1 

0.027 

731214B 

BNE 

49  8  6 

0.025 

750427B 

BNE 

18  1  1 

0.044 

760704B 

BTZ 

38  0  5 

0.030 

761123B 

BNE 

22  0  0 

0.042 

761207B 

BSW 

17  2  1 

0.044 

770329D 

Deg 

25  14  0 

0.031 

770730 D 

Deg 

21  16  0 

0.032 

780326D 

Deg 

25  6  0 

0.035 

780422D 

Deg 

21  9  0 

0.036 

4.863 


4.448 


4.791 


4.875 


5.080 


4.811 


4.481 


4.447 


4.688 


4.148 


5.619 


5.075 


5.743 


5.248 


4.904 


5.199 


4.928 


4.401 


4.296 


4.995 


4.562 


5.162 


4.766 


5.026 


5.208 


5.425 


5.073 


4.826 


4.735 


5.033 


4.405 


5.935 


5.402 


5.801 


5.985 


5.549 


5.238 


5.545 


5.351 


4.785 


4.692 


5.301 


4.821 


5.392 


5.036 


5.226 


5.442 


5.610 


5.306 


5.051 


4.991 


5.258 


4.682 


6.158 


5.624 


5.998 


6.171 


5.760 


5.521 


5.812 


5.680 


5.581 


5.073 


4.943 


5.530 


5.071 


<20 


<20 


12 


23 


20-150 


16 


6 


8 


<20 


2 


165 


20-150 


140 
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*:  a  historical  event  discussed  at  U.S.-U.S.S.R.  negotiation. 
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Event 


Magnitudes 


#  of  Signals 


Ns  Nn  Nc 


of  Semipalatinsk  Explosions 


Magnitudes  (Jih  and  Wagner, 1992) 


S.E.M. 


811129B  BSW 


811227B  BSW 


820425B  BTZ 


820704B 


820831 B  BSW 


821205B  BSW 


821226B  BNE 


83061 2B  BTZ 


831006B  BSW 


831026B  BTZ 


831120B  BNE 


840425B  BTZ 


840526 B  BNE 


840714B  BTZ 


841027B  BSW 


841 202 B  BNE 


841216B  BTZ 


841228B  BSW 


85021  OB  BSW 


85061 5B  BSW 


850630B  BSW 
870620B  BSW 
8809 14B  BSW 


37  11  5 


23  0  0 


1400 


21  0  0 


27  17  1 


26  0  0 


38  10  1 


16  0  0 


25  0  0 


1800 


1783 


21  0  0 


31  0  3 


22  0  0 


1900 


22  0  0 


15  0  0 


1900 


18  1  4 


15  0  0 


37  3  6 

24  3  13 

25  0  1 


0.027 


0.043 


0.029 


0.038 


0.028 


0.049 


0.039 


0.046 


0.037 


0.043 


0.034 


0.029 

0.031 

0.038 


4.559 


4.933 


5.547 


5.313 


5.406 

5.520 

5.480 


4.865 


5.130 


5.848 


5.679 

5.766 

5.777 


5.527 


6.196 


5.970 


6.054 


5.097 


6.108 


5.606 


5.940 


5.939 


5.989 


5.339 


5.895 


6.005 


6.057 


6.233 


5.709 


6.038 


5.924 


5.834 


6.060 


5.898 

5.999 

6.021 


:  Joint  Verification  Experiment. 


Yield  Estimates  of 


Even 

Epicenter 

Date 

Site 

Lon 

Lat 

650115B 

BTZ 

79.009 

49.935 

6511 21 D 

Deg 

78.064 

49.819 

660213D 

Deg 

78.121 

49.809 

660320D 

Deg 

78.024 

49.762 

660507D 

Deg 

78.105 

49.743 

661019D 

Deg 

78.021 

49.747 

661218M 

Mzk 

77.747 

49.925 

670226D 

Deg 

78.082 

49.746 

67091 6M 

Mzk 

77.728 

49.937 

670922M 

Mzk 

77.691 

49.960 

671122M 

Mzk 

77.687 

49.942 

68061 9B 

BNE 

78.986 

49.980 

680929D 

Deg 

78.122 

49.812 

690531 M 

Mzk 

77.694 

49.950 

690723D 

Deg 

78.130 

49.816 

69091  ID 

Deg 

77.997 

49.776 

691130B 

BTZ 

78.956 

49.924 

691228M 

Mzk 

77.714 

49.937 

700721 M 

Mzk 

77.673 

49.952 

701104M 

Mzk 

77.762 

49.989 

710322D 

Deg 

78.109 

49.798 

710425D 

Deg 

78.034 

49.769 

710606M 

Mzk 

77.660 

49.975 

Semipalatinsk  Explosions 


Yield  (Jih  and  Wagner, 


1992) 


Yield 


Announced 


94 

81 

24 

24 

197 

192 

91 

94 

1 

1 

41 

36 

75 

78 

89 

94 

10 

10 

8 

9 

1 

1 

12 

13 

53 

51 

9 

9 

13 

13 

2 

3 

101 

94 

60 

58 

13 

14 

18 

17 

40 

35 

91 

88 

23 

22 

100-150 


29 


125 


100 


20-150 


20-150 


20-150 


<20 


10 


<20 


<20 


60 


<20 


16 


<20 


125 


46 


<20 


<20 


20-150 


90 


16 
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Event 


Yield  Estimates  of 


Epicenter 


Semipalatinsk  Explosions 


Yield  (Jih  and  Wagner, 


1992) 


Date 

Site 

Lon 

Lat 

710619M 

Mzk 

77.641 

49.969 

710630B 

BTZ 

78.980 

49.946 

711009M 

Mzk 

77.641 

49.978 

711 021 M 

Mzk 

77.597 

49.974 

711230D 

Deg 

78.037 

49.760 

72021  OB 

BNE 

78.878 

50.024 

720328D 

Deg 

78.076 

49.733 

72081 6D 

Deg 

78.059 

49.765 

720826 M 

Mzk 

77.717 

49.982 

720902M 

Mzk 

77.641 

49.959 

721102B 

BSW 

78.817 

49.927 

721210D 

Deg 

78.058 

49.819 

721 21 OB 

BNE 

78.996 

50.027 

730723B 

BTZ 

78.850 

49.980 

731214B 

BNE 

79.010 

50.040 

750427B 

BNE 

78.980 

49.990 

760704B 

BTZ 

78.950 

49.910 

761123B 

BNE 

79.000 

49.990 

761207B 

BSW 

78.900 

49.880 

770329D 

Deg 

78.140 

49.790 

770730D 

Deg 

78.160 

49.770 

780326D 

Deg 

78.070 

49.730 

Yield 


Announced 


<20 


<20 


12 


23 


20-150 


16 


6 


8 


<20 


2 


165 


20-150 


140 
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Event 


Yield  Estimates  of  Semipalatinsk  Explosions 


Epicenter  Yield  (Jih  and  Wagner, 


1992) 


Date 

Site 

Lon 

Lat 

78061  IB 

BSW 

78.838 

49.879 

780705B 

BSW 

78.871 

49.887 

780728D 

Deg 

78.140 

49.756 

780829B 

BNE 

78.990 

50.000 

78091 5B 

BT2 

78.940 

49.910 

781104B 

BNE 

78.943 

50.034 

781129B 

BSW 

78.760 

49.950 

790623B 

BT2 

78.910 

49.910 

790707B 

BNE 

79.060 

50.050 

790804B 

BSW 

78.904 

49.894 

7908 18B 

BNE 

79.010 

49.970 

791028B 

BNE 

78.997 

49.973 

791202B 

BSW 

78.840 

49.890 

791223B 

BSW 

78.755 

49.916 

800522D 

Deg 

78.082 

49.784 

800629B 

BSW 

78.815 

49.939 

80091 4B 

BT2 

78.880 

49.970 

801012B 

BNE 

79.080 

49.950 

801214B 

BT2 

79.000 

49.930 

801227B 

BNE 

79.040 

50.040 

810422B 

BSW 

78.900 

49.900 

810913B 

BT2 

78.980 

49.890 

811018B 

BSW 

78.859 

49.923 

Yield 


Announced 


Event 

Epicenter  j 

Date 

Site 

Lon 

Lat 

811129B 

BSW 

78.860 

49.887 

811227B 

BSW 

78.870 

49.900 

820425B 

BTZ 

78.930 

49.880 

820704B 

BTZ 

78.850 

49.990 

820831 B 

BSW 

78.761 

49.924 

821205B 

BSW 

78.840 

49.910 

821226B 

BNE 

78.988 

50.071 

83061 2B 

BTZ 

78.980 

49.910 

831006B 

BSW 

78.840 

49.930 

831026B 

BTZ 

78.910 

49.920 

831120B 

BNE 

78.999 

50.047 

840425B 

BTZ 

78.940 

49.950 

840526 B 

BNE 

79.006 

49.969 

84071 4B 

BTZ 

78.960 

49.890 

841 027 B 

BSW 

78.830 

49.950 

841202B 

BNE 

79.070 

49.990 

841216B 

BTZ 

78.860 

49.960 

841228B 

BSW 

78.750 

49.860 

85021 OB 

BSW 

78.781 

49.888 

850615B 

BSW 

78.880 

49.890 

850630B 

BSW 

78.658 

49.848 

870620B 

BSW 

78.740 

49.927 

Yield  Estimates  of  Semipalatinsk  Explosions 


Yield  (Jih  and  Wagner,  1 992) 


Yield 


Announced 


Prof.  Thomas  Ahrens 

Seismological  Lab,  252-21 

Division  of  Geological  &  Planetary  Sciences 

California  Institute  of  Technology 

Pasadena,  C A  91125 

Prof.  Kei  iti  Aki 

Center  for  Earth  Sciences 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089-0741 

Prof.  Shelton  Alexander 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dr.  Ralph  Alewine,  HI 
DARPA/NMRO 
3701  North  Fairfax  Drive 
Arlington,  VA  22203-1714 


Prof.  Charles  B.  Archambeau 
CIRES 

University  of  Colorado 
Boulder,  CO  80309 


Dr.  Thomas  C.  Bache,  Jr. 
Science  Applicauons  Int'l  Corp. 
10260  Campus  Point  Drive 
San  Diego,  CA  92121  (2  copies) 


Prof.  Muawia  Barazangi 
Institute  for  the  Study  of  die  Continent 
Cornell  University 
Ithaca,  NY  14853 


Dr.  Jeff  Barker 

Department  of  Geological  Sciences 
State  University  of  New  York 
at  Bingharr  ton 
Vestal,  NY  13901 

Dr.  Douglas  R.  Baumgardt 
ENSCO,  Inc 
5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


Dr.  Susan  Beck 
Department  of  Geosciences 
Building  #77 
University  of  Arizona 
Tuscon,  A Z  85721 


Dr.  T.J.  Bennett 
S-CUBED 

A  Division  of  Maxwell  Laboratories 
1 1800  Sunrise  Valley  Drive,  Suite  1212 
Reston,  VA  22091 

Dr.  Robert  Blandford 
AFTAC/TT,  Center  for  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington,  VA  22209-2308 

Dr.  G.A.  Bollinger 
Department  of  Geological  Sciences 
Virginia  Polytechnical  Institute 
21044  Derring  Hall 
Blacksburg,  VA  24061 

Dr.  Stephen  Bratt 
Center  for  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington,  VA  22209-2308 

Dr.  Lawrence  Burdick 
Woodward-Clyde  Consultants 
566  El  Dorado  Street 
Pasadena,  CA  91109-3245 


Dr.  Robert  Burridge 
Schlumberger-Doll  Research  Center 
Old  Quarry  Road 
Ridgefield,  CT  06877 


Dr.  Jerry  Carter 
Center  for  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington,  VA  22209-2308 

Dr.  Eric  Chael 
Division  9241 
Sandia  Laboratory 
Albuquerque,  NM  87185 


Prof.  Vernon  F.  Cormier 
Department  of  Geology  &  Geophysics 
U-45,  Room  207 
University  of  Connecticut 
Storrs,  CT  06268 

Prof.  Steven  Day 

Department  of  Geological  Sciences 
San  Diego  State  University 
San  Diego,  CA  92182 
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Marvin  Denny 
U.S.  Department  of  Energy 
Office  of  Arms  Control 
Washington,  DC  20585 


Dr.  Zoltan  Der 
ENSCO,  Inc. 

5400  Pott  Royal  Road 
Springfield,  VA  22151-2388 


Prof.  Adam  Dziewonski 
Hoffman  Laboratory,  Harvard  University 
Dept,  of  Earth  Atmos.  &  Planetary  Sciences 
20  Oxford  Street 
Cambridge,  MA  02138 

Prof.  John  Ebel 

Department  of  Geology  &  Geophysics 
Boston  College 
Chestnut  Hill,  MA  02167 


Eric  Fielding 

SNEEHall 

INSTOC 

Cornell  University 
Ithaca,  NY  14853 

Dr.  Mark  D.  Fisk 

Mission  Research  Corporation 

735  State  Street 

P.O.  Drawer  719 

Santa  Barbara,  CA  93102 

Prof  Stanley  Flatte 
Applied  Sciences  Building 
University  of  California,  Santa  Cruz 
Santa  Cruz,  CA  95064 


Dr.  John  Foley 
NER-Geo  Sciences 
1 100  Crown  Colony  Drive 
Quincy,  MA  02169 


Prof.  Donald  Forsyth 
Department  of  Geological  Sciences 
Brown  University 
Providence,  R1  02912 


Dr.  Art  Frankel 
U.S.  Geological  Survey 
922  National  Center 
Reston,  VA  22092 


Dr.  Cliff  Frolich 
Institute  of  Geophysics 
8701  North  Mopac 
Austin,  TX  78759 


Dr.  Holly  Given 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Dr.  Jeffrey  W.  Given 
SAIC 

10260  Campus  Point  Drive 
San  Diego,  CA  92121 


Dr.  Dale  Glover 
Defense  Intelligence  Agency 
ATTN:  ODT-1B 
Washington,  DC  20301 


Dr.  Indra  Gupta 
Teledyne  Geotech 
314  Montgomery  Street 
Alexanders,  VA  22314 


Dan  N.  Hagedon 
Pacific  Northwest  Laboratories 
Battelle  Boulevard 
Richland,  WA  99352 


Dr.  James  Hannon 

Lawrence  Livermore  National  Laboratory 

P.O.  Box  808 

L-205 

Livermore,  CA  94550 

Dr.  Roger  Hansen 

HQ  AFTAQTTR 

Patrick  AFB,  FL  32925-6001 


Prof.  David  G.  Harkrider 
Seismological  Laboratory 
Division  of  Geological  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Prof.  Danny  Harvey 
CIRES 

University  of  Colorado 
Boulder,  CO  80309 
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Prof.  Donald  V.  Helmberger 
Seismological  Laboratory 
Division  of  Geological  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  C A  91125 

Prof.  Eugene  Herrin 

Institute  for  the  Study  of  Earth  and  Man 

Geophysical  Laboratory 

Southern  Methodist  University 

Dallas,  TX  75275 

Prof.  Robert  B.  Herrmann 

Department  of  Earth  &  Atmospheric  Sciences 

St.  Louis  University 

St.  Louis,  MO  63156 


Prof.  Lane  R.  Johnson 
Seismographic  Station 
University  of  California 
Berkeley,  CA  94720 


Prof.  Thomas  H.  Jordan 
Department  of  Earth,  Atmospheric  & 
Planetary  Sciences 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

Prof.  Alan  Kafka 

Department  of  Geology  &  Geophysics 
Boston  College 
Chestnut  Hill,  MA  02167 


Robert  C.  Kemerait 
ENSCO,  Inc. 

445  Pineda  Court 
Melbourne,  FL  32940 


Dr.  Max  Koontz 
U.S.  Dept,  of  Energy/DP  5 
Forrestal  Building 
1000  Independence  Avenue 
Washington,  DC  20585 

Dr.  Richard  LaCuss 

MIT  Lincoln  Laboratory,  M-200B 

P.O.  Box  73 

Lexington,  MA  02173-0073 


Dr.  Fred  K.  Lamb 

University  of  Illinois  at  Urbana-Champaign 
Department  of  Physics 
1 1 10  West  Green  Street 
Urbana,  IL  61801 


Prof.  Charles  A.  Langston 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Jim  Lawson,  Chief  Geophysicist 
Oklahoma  Geological  Survey 
Oklahoma  Geophysical  Observatory 
P.O.  Box  8 

Leonard,  OK  74043-0008 

Prof.  Thome  Lay 

Institute  of  Tectonics 

Earth  Science  Board 

University  of  California,  Santa  Cruz 

Santa  Cruz,  CA  95064 

Dr.  William  Leith 
U.S.  Geological  Survey 
Mail  Stop  928 
Reston,  VA  22092 


Mr.  James  F.  Lewkowicz 

Phillips  Laboratory/GPEH 

Hanscom  AFB,  MA  01731-5000(  2  copies) 


Mr.  Alfred  Lieberman 

ACDA/VI-OA  State  Department  Building 

Room  5726 

320-2 1st  Street,  NW 

Washington,  DC  20451 

Prof.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 


Dr.  Randolph  Martin,  HI 
New  England  Research,  Inc. 

76  Olcott  Drive 

White  River  Junction,  VT  05001 


Dr.  Robert  Masse 
Denver  Federal  Building 
Box  25046,  Mail  Stop  967 
Denver,  CO  80225 


Dr.  Gary  McCartor 
Department  of  Physics 
Southern  Methodist  University 
Dallas,  TX  75275 
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